\(\int \frac {\cos ^4(x)}{a+b \cos (x)+c \cos ^2(x)} \, dx\) [13]

Optimal result

Integrand size = 19, antiderivative size = 326 \[ \int \frac {\cos ^4(x)}{a+b \cos (x)+c \cos ^2(x)} \, dx=\frac {x}{2 c}+\frac {\left (b^2-a c\right ) x}{c^3}-\frac {2 \left (b^3-2 a b c-\frac {b^4-4 a b^2 c+2 a^2 c^2}{\sqrt {b^2-4 a c}}\right ) \arctan \left (\frac {\sqrt {b-2 c-\sqrt {b^2-4 a c}} \tan \left (\frac {x}{2}\right )}{\sqrt {b+2 c-\sqrt {b^2-4 a c}}}\right )}{c^3 \sqrt {b-2 c-\sqrt {b^2-4 a c}} \sqrt {b+2 c-\sqrt {b^2-4 a c}}}-\frac {2 \left (b^3-2 a b c+\frac {b^4-4 a b^2 c+2 a^2 c^2}{\sqrt {b^2-4 a c}}\right ) \arctan \left (\frac {\sqrt {b-2 c+\sqrt {b^2-4 a c}} \tan \left (\frac {x}{2}\right )}{\sqrt {b+2 c+\sqrt {b^2-4 a c}}}\right )}{c^3 \sqrt {b-2 c+\sqrt {b^2-4 a c}} \sqrt {b+2 c+\sqrt {b^2-4 a c}}}-\frac {b \sin (x)}{c^2}+\frac {\cos (x) \sin (x)}{2 c} \]

[Out]

1/2*x/c+(-a*c+b^2)*x/c^3-b*sin(x)/c^2+1/2*cos(x)*sin(x)/c-2*arctan((b-2*c-(-4*a*c+b^2)^(1/2))^(1/2)*tan(1/2*x)
/(b+2*c-(-4*a*c+b^2)^(1/2))^(1/2))*(b^3-2*a*b*c+(-2*a^2*c^2+4*a*b^2*c-b^4)/(-4*a*c+b^2)^(1/2))/c^3/(b-2*c-(-4*
a*c+b^2)^(1/2))^(1/2)/(b+2*c-(-4*a*c+b^2)^(1/2))^(1/2)-2*arctan((b-2*c+(-4*a*c+b^2)^(1/2))^(1/2)*tan(1/2*x)/(b
+2*c+(-4*a*c+b^2)^(1/2))^(1/2))*(b^3-2*a*b*c+(2*a^2*c^2-4*a*b^2*c+b^4)/(-4*a*c+b^2)^(1/2))/c^3/(b-2*c+(-4*a*c+
b^2)^(1/2))^(1/2)/(b+2*c+(-4*a*c+b^2)^(1/2))^(1/2)

Rubi [A] (verified)

Time = 4.62 (sec) , antiderivative size = 326, normalized size of antiderivative = 1.00, number of steps used = 10, number of rules used = 7, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.368, Rules used = {3338, 2717, 2715, 8, 3374, 2738, 211} \[ \int \frac {\cos ^4(x)}{a+b \cos (x)+c \cos ^2(x)} \, dx=-\frac {2 \left (-\frac {2 a^2 c^2-4 a b^2 c+b^4}{\sqrt {b^2-4 a c}}-2 a b c+b^3\right ) \arctan \left (\frac {\tan \left (\frac {x}{2}\right ) \sqrt {-\sqrt {b^2-4 a c}+b-2 c}}{\sqrt {-\sqrt {b^2-4 a c}+b+2 c}}\right )}{c^3 \sqrt {-\sqrt {b^2-4 a c}+b-2 c} \sqrt {-\sqrt {b^2-4 a c}+b+2 c}}-\frac {2 \left (\frac {2 a^2 c^2-4 a b^2 c+b^4}{\sqrt {b^2-4 a c}}-2 a b c+b^3\right ) \arctan \left (\frac {\tan \left (\frac {x}{2}\right ) \sqrt {\sqrt {b^2-4 a c}+b-2 c}}{\sqrt {\sqrt {b^2-4 a c}+b+2 c}}\right )}{c^3 \sqrt {\sqrt {b^2-4 a c}+b-2 c} \sqrt {\sqrt {b^2-4 a c}+b+2 c}}+\frac {x \left (b^2-a c\right )}{c^3}-\frac {b \sin (x)}{c^2}+\frac {x}{2 c}+\frac {\sin (x) \cos (x)}{2 c} \]

[In]

Int[Cos[x]^4/(a + b*Cos[x] + c*Cos[x]^2),x]

[Out]

x/(2*c) + ((b^2 - a*c)*x)/c^3 - (2*(b^3 - 2*a*b*c - (b^4 - 4*a*b^2*c + 2*a^2*c^2)/Sqrt[b^2 - 4*a*c])*ArcTan[(S
qrt[b - 2*c - Sqrt[b^2 - 4*a*c]]*Tan[x/2])/Sqrt[b + 2*c - Sqrt[b^2 - 4*a*c]]])/(c^3*Sqrt[b - 2*c - Sqrt[b^2 -
4*a*c]]*Sqrt[b + 2*c - Sqrt[b^2 - 4*a*c]]) - (2*(b^3 - 2*a*b*c + (b^4 - 4*a*b^2*c + 2*a^2*c^2)/Sqrt[b^2 - 4*a*
c])*ArcTan[(Sqrt[b - 2*c + Sqrt[b^2 - 4*a*c]]*Tan[x/2])/Sqrt[b + 2*c + Sqrt[b^2 - 4*a*c]]])/(c^3*Sqrt[b - 2*c
+ Sqrt[b^2 - 4*a*c]]*Sqrt[b + 2*c + Sqrt[b^2 - 4*a*c]]) - (b*Sin[x])/c^2 + (Cos[x]*Sin[x])/(2*c)

Rule 8

Int[a_, x_Symbol] :> Simp[a*x, x] /; FreeQ[a, x]

Rule 211

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]/a)*ArcTan[x/Rt[a/b, 2]], x] /; FreeQ[{a, b}, x]
&& PosQ[a/b]

Rule 2715

Int[((b_.)*sin[(c_.) + (d_.)*(x_)])^(n_), x_Symbol] :> Simp[(-b)*Cos[c + d*x]*((b*Sin[c + d*x])^(n - 1)/(d*n))
, x] + Dist[b^2*((n - 1)/n), Int[(b*Sin[c + d*x])^(n - 2), x], x] /; FreeQ[{b, c, d}, x] && GtQ[n, 1] && Integ
erQ[2*n]

Rule 2717

Int[sin[Pi/2 + (c_.) + (d_.)*(x_)], x_Symbol] :> Simp[Sin[c + d*x]/d, x] /; FreeQ[{c, d}, x]

Rule 2738

Int[((a_) + (b_.)*sin[Pi/2 + (c_.) + (d_.)*(x_)])^(-1), x_Symbol] :> With[{e = FreeFactors[Tan[(c + d*x)/2], x
]}, Dist[2*(e/d), Subst[Int[1/(a + b + (a - b)*e^2*x^2), x], x, Tan[(c + d*x)/2]/e], x]] /; FreeQ[{a, b, c, d}
, x] && NeQ[a^2 - b^2, 0]

Rule 3338

Int[cos[(d_.) + (e_.)*(x_)]^(m_.)*((a_.) + cos[(d_.) + (e_.)*(x_)]^(n_.)*(b_.) + cos[(d_.) + (e_.)*(x_)]^(n2_.
)*(c_.))^(p_), x_Symbol] :> Int[ExpandTrig[cos[d + e*x]^m*(a + b*cos[d + e*x]^n + c*cos[d + e*x]^(2*n))^p, x],
x] /; FreeQ[{a, b, c, d, e}, x] && EqQ[n2, 2*n] && NeQ[b^2 - 4*a*c, 0] && IntegersQ[m, n, p]

Rule 3374

Int[(cos[(d_.) + (e_.)*(x_)]*(B_.) + (A_))/((a_.) + cos[(d_.) + (e_.)*(x_)]*(b_.) + cos[(d_.) + (e_.)*(x_)]^2*
(c_.)), x_Symbol] :> Module[{q = Rt[b^2 - 4*a*c, 2]}, Dist[B + (b*B - 2*A*c)/q, Int[1/(b + q + 2*c*Cos[d + e*x
]), x], x] + Dist[B - (b*B - 2*A*c)/q, Int[1/(b - q + 2*c*Cos[d + e*x]), x], x]] /; FreeQ[{a, b, c, d, e, A, B
}, x] && NeQ[b^2 - 4*a*c, 0]

Rubi steps \begin{align*} \text {integral}& = \int \left (\frac {b^2-a c}{c^3}-\frac {b \cos (x)}{c^2}+\frac {\cos ^2(x)}{c}+\frac {-a b^2 \left (1-\frac {a c}{b^2}\right )-b^3 \left (1-\frac {2 a c}{b^2}\right ) \cos (x)}{c^3 \left (a+b \cos (x)+c \cos ^2(x)\right )}\right ) \, dx \\ & = \frac {\left (b^2-a c\right ) x}{c^3}+\frac {\int \frac {-a b^2 \left (1-\frac {a c}{b^2}\right )-b^3 \left (1-\frac {2 a c}{b^2}\right ) \cos (x)}{a+b \cos (x)+c \cos ^2(x)} \, dx}{c^3}-\frac {b \int \cos (x) \, dx}{c^2}+\frac {\int \cos ^2(x) \, dx}{c} \\ & = \frac {\left (b^2-a c\right ) x}{c^3}-\frac {b \sin (x)}{c^2}+\frac {\cos (x) \sin (x)}{2 c}+\frac {\int 1 \, dx}{2 c}-\frac {\left (b^3-2 a b c-\frac {b^4-4 a b^2 c+2 a^2 c^2}{\sqrt {b^2-4 a c}}\right ) \int \frac {1}{b-\sqrt {b^2-4 a c}+2 c \cos (x)} \, dx}{c^3}-\frac {\left (b^3-2 a b c+\frac {b^4-4 a b^2 c+2 a^2 c^2}{\sqrt {b^2-4 a c}}\right ) \int \frac {1}{b+\sqrt {b^2-4 a c}+2 c \cos (x)} \, dx}{c^3} \\ & = \frac {x}{2 c}+\frac {\left (b^2-a c\right ) x}{c^3}-\frac {b \sin (x)}{c^2}+\frac {\cos (x) \sin (x)}{2 c}-\frac {\left (2 \left (b^3-2 a b c-\frac {b^4-4 a b^2 c+2 a^2 c^2}{\sqrt {b^2-4 a c}}\right )\right ) \text {Subst}\left (\int \frac {1}{b+2 c-\sqrt {b^2-4 a c}+\left (b-2 c-\sqrt {b^2-4 a c}\right ) x^2} \, dx,x,\tan \left (\frac {x}{2}\right )\right )}{c^3}-\frac {\left (2 \left (b^3-2 a b c+\frac {b^4-4 a b^2 c+2 a^2 c^2}{\sqrt {b^2-4 a c}}\right )\right ) \text {Subst}\left (\int \frac {1}{b+2 c+\sqrt {b^2-4 a c}+\left (b-2 c+\sqrt {b^2-4 a c}\right ) x^2} \, dx,x,\tan \left (\frac {x}{2}\right )\right )}{c^3} \\ & = \frac {x}{2 c}+\frac {\left (b^2-a c\right ) x}{c^3}-\frac {2 \left (b^3-2 a b c-\frac {b^4-4 a b^2 c+2 a^2 c^2}{\sqrt {b^2-4 a c}}\right ) \arctan \left (\frac {\sqrt {b-2 c-\sqrt {b^2-4 a c}} \tan \left (\frac {x}{2}\right )}{\sqrt {b+2 c-\sqrt {b^2-4 a c}}}\right )}{c^3 \sqrt {b-2 c-\sqrt {b^2-4 a c}} \sqrt {b+2 c-\sqrt {b^2-4 a c}}}-\frac {2 \left (b^3-2 a b c+\frac {b^4-4 a b^2 c+2 a^2 c^2}{\sqrt {b^2-4 a c}}\right ) \arctan \left (\frac {\sqrt {b-2 c+\sqrt {b^2-4 a c}} \tan \left (\frac {x}{2}\right )}{\sqrt {b+2 c+\sqrt {b^2-4 a c}}}\right )}{c^3 \sqrt {b-2 c+\sqrt {b^2-4 a c}} \sqrt {b+2 c+\sqrt {b^2-4 a c}}}-\frac {b \sin (x)}{c^2}+\frac {\cos (x) \sin (x)}{2 c} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.60 (sec) , antiderivative size = 356, normalized size of antiderivative = 1.09 \[ \int \frac {\cos ^4(x)}{a+b \cos (x)+c \cos ^2(x)} \, dx=\frac {4 b^2 x+2 c (-2 a+c) x+\frac {4 \sqrt {2} \left (b^4-4 a b^2 c+2 a^2 c^2+b^3 \sqrt {b^2-4 a c}-2 a b c \sqrt {b^2-4 a c}\right ) \text {arctanh}\left (\frac {\left (b-2 c+\sqrt {b^2-4 a c}\right ) \tan \left (\frac {x}{2}\right )}{\sqrt {-2 b^2+4 c (a+c)-2 b \sqrt {b^2-4 a c}}}\right )}{\sqrt {b^2-4 a c} \sqrt {-b^2+2 c (a+c)-b \sqrt {b^2-4 a c}}}-\frac {4 \sqrt {2} \left (-b^4+4 a b^2 c-2 a^2 c^2+b^3 \sqrt {b^2-4 a c}-2 a b c \sqrt {b^2-4 a c}\right ) \text {arctanh}\left (\frac {\left (-b+2 c+\sqrt {b^2-4 a c}\right ) \tan \left (\frac {x}{2}\right )}{\sqrt {-2 b^2+4 c (a+c)+2 b \sqrt {b^2-4 a c}}}\right )}{\sqrt {b^2-4 a c} \sqrt {-b^2+2 c (a+c)+b \sqrt {b^2-4 a c}}}-4 b c \sin (x)+c^2 \sin (2 x)}{4 c^3} \]

[In]

Integrate[Cos[x]^4/(a + b*Cos[x] + c*Cos[x]^2),x]

[Out]

(4*b^2*x + 2*c*(-2*a + c)*x + (4*Sqrt[2]*(b^4 - 4*a*b^2*c + 2*a^2*c^2 + b^3*Sqrt[b^2 - 4*a*c] - 2*a*b*c*Sqrt[b
^2 - 4*a*c])*ArcTanh[((b - 2*c + Sqrt[b^2 - 4*a*c])*Tan[x/2])/Sqrt[-2*b^2 + 4*c*(a + c) - 2*b*Sqrt[b^2 - 4*a*c
]]])/(Sqrt[b^2 - 4*a*c]*Sqrt[-b^2 + 2*c*(a + c) - b*Sqrt[b^2 - 4*a*c]]) - (4*Sqrt[2]*(-b^4 + 4*a*b^2*c - 2*a^2
*c^2 + b^3*Sqrt[b^2 - 4*a*c] - 2*a*b*c*Sqrt[b^2 - 4*a*c])*ArcTanh[((-b + 2*c + Sqrt[b^2 - 4*a*c])*Tan[x/2])/Sq
rt[-2*b^2 + 4*c*(a + c) + 2*b*Sqrt[b^2 - 4*a*c]]])/(Sqrt[b^2 - 4*a*c]*Sqrt[-b^2 + 2*c*(a + c) + b*Sqrt[b^2 - 4
*a*c]]) - 4*b*c*Sin[x] + c^2*Sin[2*x])/(4*c^3)

Maple [A] (verified)

Time = 4.72 (sec) , antiderivative size = 431, normalized size of antiderivative = 1.32

method	result	size
default	\(\frac {2 \left (a -b +c \right ) \left (\frac {\left (a^{2} c \sqrt {-4 a c +b^{2}}-a \,b^{2} \sqrt {-4 a c +b^{2}}-2 c a b \sqrt {-4 a c +b^{2}}+b^{3} \sqrt {-4 a c +b^{2}}+3 c \,a^{2} b +2 a^{2} c^{2}-a \,b^{3}-4 a \,b^{2} c +b^{4}\right ) \operatorname {arctanh}\left (\frac {\left (-a +b -c \right ) \tan \left (\frac {x}{2}\right )}{\sqrt {\left (\sqrt {-4 a c +b^{2}}-a +c \right ) \left (a -b +c \right )}}\right )}{2 \sqrt {-4 a c +b^{2}}\, \left (a -b +c \right ) \sqrt {\left (\sqrt {-4 a c +b^{2}}-a +c \right ) \left (a -b +c \right )}}+\frac {\left (a^{2} c \sqrt {-4 a c +b^{2}}-a \,b^{2} \sqrt {-4 a c +b^{2}}-2 c a b \sqrt {-4 a c +b^{2}}+b^{3} \sqrt {-4 a c +b^{2}}-3 c \,a^{2} b -2 a^{2} c^{2}+a \,b^{3}+4 a \,b^{2} c -b^{4}\right ) \arctan \left (\frac {\left (a -b +c \right ) \tan \left (\frac {x}{2}\right )}{\sqrt {\left (\sqrt {-4 a c +b^{2}}+a -c \right ) \left (a -b +c \right )}}\right )}{2 \sqrt {-4 a c +b^{2}}\, \left (a -b +c \right ) \sqrt {\left (\sqrt {-4 a c +b^{2}}+a -c \right ) \left (a -b +c \right )}}\right )}{c^{3}}-\frac {2 \left (\frac {\left (c b +\frac {1}{2} c^{2}\right ) \left (\tan ^{3}\left (\frac {x}{2}\right )\right )+\left (c b -\frac {1}{2} c^{2}\right ) \tan \left (\frac {x}{2}\right )}{\left (1+\tan ^{2}\left (\frac {x}{2}\right )\right )^{2}}+\frac {\left (2 a c -2 b^{2}-c^{2}\right ) \arctan \left (\tan \left (\frac {x}{2}\right )\right )}{2}\right )}{c^{3}}\)	\(431\)
risch	\(\text {Expression too large to display}\)	\(3886\)

[In]

int(cos(x)^4/(a+cos(x)*b+c*cos(x)^2),x,method=_RETURNVERBOSE)

[Out]

2/c^3*(a-b+c)*(1/2*(a^2*c*(-4*a*c+b^2)^(1/2)-a*b^2*(-4*a*c+b^2)^(1/2)-2*c*a*b*(-4*a*c+b^2)^(1/2)+b^3*(-4*a*c+b
^2)^(1/2)+3*c*a^2*b+2*a^2*c^2-a*b^3-4*a*b^2*c+b^4)/(-4*a*c+b^2)^(1/2)/(a-b+c)/(((-4*a*c+b^2)^(1/2)-a+c)*(a-b+c
))^(1/2)*arctanh((-a+b-c)*tan(1/2*x)/(((-4*a*c+b^2)^(1/2)-a+c)*(a-b+c))^(1/2))+1/2*(a^2*c*(-4*a*c+b^2)^(1/2)-a
*b^2*(-4*a*c+b^2)^(1/2)-2*c*a*b*(-4*a*c+b^2)^(1/2)+b^3*(-4*a*c+b^2)^(1/2)-3*c*a^2*b-2*a^2*c^2+a*b^3+4*a*b^2*c-
b^4)/(-4*a*c+b^2)^(1/2)/(a-b+c)/(((-4*a*c+b^2)^(1/2)+a-c)*(a-b+c))^(1/2)*arctan((a-b+c)*tan(1/2*x)/(((-4*a*c+b
^2)^(1/2)+a-c)*(a-b+c))^(1/2)))-2/c^3*(((c*b+1/2*c^2)*tan(1/2*x)^3+(c*b-1/2*c^2)*tan(1/2*x))/(1+tan(1/2*x)^2)^
2+1/2*(2*a*c-2*b^2-c^2)*arctan(tan(1/2*x)))

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 8167 vs. \(2 (282) = 564\).

Time = 5.52 (sec) , antiderivative size = 8167, normalized size of antiderivative = 25.05 \[ \int \frac {\cos ^4(x)}{a+b \cos (x)+c \cos ^2(x)} \, dx=\text {Too large to display} \]

[In]

integrate(cos(x)^4/(a+b*cos(x)+c*cos(x)^2),x, algorithm="fricas")

[Out]

Too large to include

Sympy [F(-1)]

Timed out. \[ \int \frac {\cos ^4(x)}{a+b \cos (x)+c \cos ^2(x)} \, dx=\text {Timed out} \]

[In]

integrate(cos(x)**4/(a+b*cos(x)+c*cos(x)**2),x)

[Out]

Timed out

Maxima [F]

\[ \int \frac {\cos ^4(x)}{a+b \cos (x)+c \cos ^2(x)} \, dx=\int { \frac {\cos \left (x\right )^{4}}{c \cos \left (x\right )^{2} + b \cos \left (x\right ) + a} \,d x } \]

[In]

integrate(cos(x)^4/(a+b*cos(x)+c*cos(x)^2),x, algorithm="maxima")

[Out]

1/4*(4*c^3*integrate(-2*(2*(b^4 - 2*a*b^2*c)*cos(3*x)^2 + 4*(2*a^2*b^2 - a^2*c^2 - (2*a^3 - a*b^2)*c)*cos(2*x)
^2 + 2*(b^4 - 2*a*b^2*c)*cos(x)^2 + 2*(b^4 - 2*a*b^2*c)*sin(3*x)^2 + 4*(2*a^2*b^2 - a^2*c^2 - (2*a^3 - a*b^2)*
c)*sin(2*x)^2 + 2*(4*a*b^3 - 2*a*b*c^2 - (6*a^2*b - b^3)*c)*sin(2*x)*sin(x) + 2*(b^4 - 2*a*b^2*c)*sin(x)^2 + (
(b^3*c - 2*a*b*c^2)*cos(3*x) + 2*(a*b^2*c - a^2*c^2)*cos(2*x) + (b^3*c - 2*a*b*c^2)*cos(x))*cos(4*x) + (b^3*c
- 2*a*b*c^2 + 2*(4*a*b^3 - 2*a*b*c^2 - (6*a^2*b - b^3)*c)*cos(2*x) + 4*(b^4 - 2*a*b^2*c)*cos(x))*cos(3*x) + 2*
(a*b^2*c - a^2*c^2 + (4*a*b^3 - 2*a*b*c^2 - (6*a^2*b - b^3)*c)*cos(x))*cos(2*x) + (b^3*c - 2*a*b*c^2)*cos(x) +
((b^3*c - 2*a*b*c^2)*sin(3*x) + 2*(a*b^2*c - a^2*c^2)*sin(2*x) + (b^3*c - 2*a*b*c^2)*sin(x))*sin(4*x) + 2*((4
*a*b^3 - 2*a*b*c^2 - (6*a^2*b - b^3)*c)*sin(2*x) + 2*(b^4 - 2*a*b^2*c)*sin(x))*sin(3*x))/(c^5*cos(4*x)^2 + 4*b
^2*c^3*cos(3*x)^2 + 4*b^2*c^3*cos(x)^2 + c^5*sin(4*x)^2 + 4*b^2*c^3*sin(3*x)^2 + 4*b^2*c^3*sin(x)^2 + 4*b*c^4*
cos(x) + c^5 + 4*(4*a^2*c^3 + 4*a*c^4 + c^5)*cos(2*x)^2 + 4*(4*a^2*c^3 + 4*a*c^4 + c^5)*sin(2*x)^2 + 8*(2*a*b*
c^3 + b*c^4)*sin(2*x)*sin(x) + 2*(2*b*c^4*cos(3*x) + 2*b*c^4*cos(x) + c^5 + 2*(2*a*c^4 + c^5)*cos(2*x))*cos(4*
x) + 4*(2*b^2*c^3*cos(x) + b*c^4 + 2*(2*a*b*c^3 + b*c^4)*cos(2*x))*cos(3*x) + 4*(2*a*c^4 + c^5 + 2*(2*a*b*c^3
+ b*c^4)*cos(x))*cos(2*x) + 4*(b*c^4*sin(3*x) + b*c^4*sin(x) + (2*a*c^4 + c^5)*sin(2*x))*sin(4*x) + 8*(b^2*c^3
*sin(x) + (2*a*b*c^3 + b*c^4)*sin(2*x))*sin(3*x)), x) + c^2*sin(2*x) - 4*b*c*sin(x) + 2*(2*b^2 - 2*a*c + c^2)*
x)/c^3

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 12587 vs. \(2 (282) = 564\).

Time = 2.78 (sec) , antiderivative size = 12587, normalized size of antiderivative = 38.61 \[ \int \frac {\cos ^4(x)}{a+b \cos (x)+c \cos ^2(x)} \, dx=\text {Too large to display} \]

[In]

integrate(cos(x)^4/(a+b*cos(x)+c*cos(x)^2),x, algorithm="giac")

[Out]

-((2*a^3*b^6 - 6*a^2*b^7 + 6*a*b^8 - 2*b^9 - 18*a^4*b^4*c + 56*a^3*b^5*c - 54*a^2*b^6*c + 12*a*b^7*c + 4*b^8*c
+ 48*a^5*b^2*c^2 - 160*a^4*b^3*c^2 + 140*a^3*b^4*c^2 + 12*a^2*b^5*c^2 - 38*a*b^6*c^2 - 2*b^7*c^2 - 32*a^6*c^3
+ 128*a^5*b*c^3 - 64*a^4*b^2*c^3 - 160*a^3*b^3*c^3 + 110*a^2*b^4*c^3 + 20*a*b^5*c^3 - 64*a^5*c^4 + 192*a^4*b*
c^4 - 80*a^3*b^2*c^4 - 64*a^2*b^3*c^4 - 32*a^4*c^5 + 64*a^3*b*c^5 + 3*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 -
4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^3*b^4 - 2*(b^2 - 4*a*c)*a^3*b^4 - 5*sqrt(a^2 - a*b + b*c - c^2 + sqrt(
b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^2*b^5 + 6*(b^2 - 4*a*c)*a^2*b^5 - 3*sqrt(a^2 - a*b + b*c - c^2 +
sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a*b^6 - 6*(b^2 - 4*a*c)*a*b^6 + 5*sqrt(a^2 - a*b + b*c - c^2
+ sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*b^7 + 2*(b^2 - 4*a*c)*b^7 - 15*sqrt(a^2 - a*b + b*c - c^2
+ sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^4*b^2*c + 10*(b^2 - 4*a*c)*a^4*b^2*c + 28*sqrt(a^2 - a*b
+ b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^3*b^3*c - 32*(b^2 - 4*a*c)*a^3*b^3*c + 27*sqr
t(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^2*b^4*c + 30*(b^2 - 4*a*c)*a^2*b^
4*c - 38*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a*b^5*c - 4*(b^2 - 4*a*
c)*a*b^5*c - 6*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*b^6*c - 4*(b^2 -
4*a*c)*b^6*c + 12*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^5*c^2 - 8*(b
^2 - 4*a*c)*a^5*c^2 - 32*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^4*b*c
^2 + 32*(b^2 - 4*a*c)*a^4*b*c^2 - 74*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*
a*c)*a^3*b^2*c^2 - 20*(b^2 - 4*a*c)*a^3*b^2*c^2 + 94*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c
))*sqrt(b^2 - 4*a*c)*a^2*b^3*c^2 - 28*(b^2 - 4*a*c)*a^2*b^3*c^2 + 31*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4
*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a*b^4*c^2 + 22*(b^2 - 4*a*c)*a*b^4*c^2 + 5*sqrt(a^2 - a*b + b*c - c^2 + s
qrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*b^5*c^2 + 2*(b^2 - 4*a*c)*b^5*c^2 + 56*sqrt(a^2 - a*b + b*c -
c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^4*c^3 - 16*(b^2 - 4*a*c)*a^4*c^3 - 88*sqrt(a^2 - a*b
+ b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^3*b*c^3 + 48*(b^2 - 4*a*c)*a^3*b*c^3 - 23*sqr
t(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^2*b^2*c^3 - 22*(b^2 - 4*a*c)*a^2*
b^2*c^3 - 30*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a*b^3*c^3 - 12*(b^2
- 4*a*c)*a*b^3*c^3 - 20*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^3*c^4
- 8*(b^2 - 4*a*c)*a^3*c^4 + 40*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*
a^2*b*c^4 + 16*(b^2 - 4*a*c)*a^2*b*c^4)*c^2*abs(a - b + c) + (4*a^3*b^6*c - 4*a^2*b^7*c - 4*a*b^8*c + 4*b^9*c
- 36*a^4*b^4*c^2 + 32*a^3*b^5*c^2 + 52*a^2*b^6*c^2 - 40*a*b^7*c^2 - 8*b^8*c^2 + 96*a^5*b^2*c^3 - 64*a^4*b^3*c^
3 - 232*a^3*b^4*c^3 + 120*a^2*b^5*c^3 + 84*a*b^6*c^3 + 4*b^7*c^3 - 64*a^6*c^4 + 384*a^4*b^2*c^4 - 64*a^3*b^3*c
^4 - 292*a^2*b^4*c^4 - 40*a*b^5*c^4 - 128*a^5*c^5 - 128*a^4*b*c^5 + 352*a^3*b^2*c^5 + 128*a^2*b^3*c^5 - 64*a^4
*c^6 - 128*a^3*b*c^6 + 3*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*a^4*b^4*c - 2*sqrt(a^2 -
a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*a^3*b^5*c - 8*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*
(a - b + c))*a^2*b^6*c + 2*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*a*b^7*c + 5*sqrt(a^2 -
a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*b^8*c - 15*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a
- b + c))*a^5*b^2*c^2 + 7*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*a^4*b^3*c^2 + 62*sqrt(a^
2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*a^3*b^4*c^2 - 51*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 -
4*a*c)*(a - b + c))*a*b^6*c^2 - 11*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*b^7*c^2 + 12*sq
rt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*a^6*c^3 + 4*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 -
4*a*c)*(a - b + c))*a^5*b*c^3 - 137*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*a^4*b^2*c^3 -
56*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*a^3*b^3*c^3 + 164*sqrt(a^2 - a*b + b*c - c^2 +
sqrt(b^2 - 4*a*c)*(a - b + c))*a^2*b^4*c^3 + 86*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*a*
b^5*c^3 + 11*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*b^6*c^3 + 68*sqrt(a^2 - a*b + b*c - c
^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*a^5*c^4 + 96*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*a
^4*b*c^4 - 169*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*a^3*b^2*c^4 - 197*sqrt(a^2 - a*b +
b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*a^2*b^3*c^4 - 71*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a
- b + c))*a*b^4*c^4 - 5*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*b^5*c^4 + 36*sqrt(a^2 - a*
b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*a^4*c^5 + 116*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a
- b + c))*a^3*b*c^5 + 113*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*a^2*b^2*c^5 + 30*sqrt(a
^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*a*b^3*c^5 - 20*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4
*a*c)*(a - b + c))*a^3*c^6 - 40*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*a^2*b*c^6 - 4*(b^2
- 4*a*c)*a^3*b^4*c + 4*(b^2 - 4*a*c)*a^2*b^5*c + 4*(b^2 - 4*a*c)*a*b^6*c - 4*(b^2 - 4*a*c)*b^7*c + 20*(b^2 -
4*a*c)*a^4*b^2*c^2 - 16*(b^2 - 4*a*c)*a^3*b^3*c^2 - 36*(b^2 - 4*a*c)*a^2*b^4*c^2 + 24*(b^2 - 4*a*c)*a*b^5*c^2
+ 8*(b^2 - 4*a*c)*b^6*c^2 - 16*(b^2 - 4*a*c)*a^5*c^3 + 88*(b^2 - 4*a*c)*a^3*b^2*c^3 - 24*(b^2 - 4*a*c)*a^2*b^3
*c^3 - 52*(b^2 - 4*a*c)*a*b^4*c^3 - 4*(b^2 - 4*a*c)*b^5*c^3 - 32*(b^2 - 4*a*c)*a^4*c^4 - 32*(b^2 - 4*a*c)*a^3*
b*c^4 + 84*(b^2 - 4*a*c)*a^2*b^2*c^4 + 24*(b^2 - 4*a*c)*a*b^3*c^4 - 16*(b^2 - 4*a*c)*a^3*c^5 - 32*(b^2 - 4*a*c
)*a^2*b*c^5)*abs(a - b + c)*abs(c) + (2*a^4*b^5*c^2 - 6*a^3*b^6*c^2 + 6*a^2*b^7*c^2 - 2*a*b^8*c^2 - 14*a^5*b^3
*c^3 + 44*a^4*b^4*c^3 - 44*a^3*b^5*c^3 + 14*a^2*b^6*c^3 - 2*a*b^7*c^3 + 2*b^8*c^3 + 24*a^6*b*c^4 - 84*a^5*b^2*
c^4 + 82*a^4*b^3*c^4 - 20*a^3*b^4*c^4 + 12*a^2*b^5*c^4 - 10*a*b^6*c^4 - 4*b^7*c^4 + 16*a^6*c^5 - 8*a^5*b*c^5 -
20*a^4*b^2*c^5 - 10*a^3*b^3*c^5 - 8*a^2*b^4*c^5 + 30*a*b^5*c^5 + 2*b^6*c^5 + 16*a^5*c^6 - 24*a^4*b*c^6 + 68*a
^3*b^2*c^6 - 58*a^2*b^3*c^6 - 16*a*b^4*c^6 - 16*a^4*c^7 + 8*a^3*b*c^7 + 36*a^2*b^2*c^7 - 16*a^3*c^8 + 3*sqrt(a
^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^4*b^3*c^2 - 2*(b^2 - 4*a*c)*a^4*b^3*
c^2 - 5*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^3*b^4*c^2 + 6*(b^2 - 4
*a*c)*a^3*b^4*c^2 - 3*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^2*b^5*c^
2 - 6*(b^2 - 4*a*c)*a^2*b^5*c^2 + 5*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a
*c)*a*b^6*c^2 + 2*(b^2 - 4*a*c)*a*b^6*c^2 - 9*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt
(b^2 - 4*a*c)*a^5*b*c^3 + 6*(b^2 - 4*a*c)*a^5*b*c^3 + 18*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b
+ c))*sqrt(b^2 - 4*a*c)*a^4*b^2*c^3 - 20*(b^2 - 4*a*c)*a^4*b^2*c^3 + 18*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2
- 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^3*b^3*c^3 + 20*(b^2 - 4*a*c)*a^3*b^3*c^3 - 23*sqrt(a^2 - a*b + b*c
- c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^2*b^4*c^3 - 6*(b^2 - 4*a*c)*a^2*b^4*c^3 - 3*sqrt(a^
2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a*b^5*c^3 + 2*(b^2 - 4*a*c)*a*b^5*c^3 -
5*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*b^6*c^3 - 2*(b^2 - 4*a*c)*b^6
*c^3 - 6*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^5*c^4 + 4*(b^2 - 4*a*
c)*a^5*c^4 - 29*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^4*b*c^4 - 2*(b
^2 - 4*a*c)*a^4*b*c^4 + 30*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^3*b
^2*c^4 - 4*(b^2 - 4*a*c)*a^3*b^2*c^4 - 2*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2
- 4*a*c)*a^2*b^3*c^4 - 4*(b^2 - 4*a*c)*a^2*b^3*c^4 + 33*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b
+ c))*sqrt(b^2 - 4*a*c)*a*b^4*c^4 + 2*(b^2 - 4*a*c)*a*b^4*c^4 + 6*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*
c)*(a - b + c))*sqrt(b^2 - 4*a*c)*b^5*c^4 + 4*(b^2 - 4*a*c)*b^5*c^4 - 22*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2
- 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^4*c^5 + 4*(b^2 - 4*a*c)*a^4*c^5 + 41*sqrt(a^2 - a*b + b*c - c^2 + s
qrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^3*b*c^5 - 6*(b^2 - 4*a*c)*a^3*b*c^5 - 68*sqrt(a^2 - a*b + b*
c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^2*b^2*c^5 + 16*(b^2 - 4*a*c)*a^2*b^2*c^5 - 19*sqr
t(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a*b^3*c^5 - 14*(b^2 - 4*a*c)*a*b^3*
c^5 - 5*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*b^4*c^5 - 2*(b^2 - 4*a*c
)*b^4*c^5 + 38*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^3*c^6 - 4*(b^2
- 4*a*c)*a^3*c^6 - 3*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^2*b*c^6 +
2*(b^2 - 4*a*c)*a^2*b*c^6 + 20*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*
a*b^2*c^6 + 8*(b^2 - 4*a*c)*a*b^2*c^6 - 10*sqrt(a^2 - a*b + b*c - c^2 + sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^
2 - 4*a*c)*a^2*c^7 - 4*(b^2 - 4*a*c)*a^2*c^7)*abs(a - b + c))*(pi*floor(1/2*x/pi + 1/2) + arctan(2*sqrt(1/2)*t
an(1/2*x)/sqrt((2*a*c^3 - 2*c^4 + sqrt(-4*(a*c^3 + b*c^3 + c^4)*(a*c^3 - b*c^3 + c^4) + 4*(a*c^3 - c^4)^2))/(a
*c^3 - b*c^3 + c^4))))/((3*a^6*b^2*c^4 - 8*a^5*b^3*c^4 - a^4*b^4*c^4 + 16*a^3*b^5*c^4 - 7*a^2*b^6*c^4 - 8*a*b^
7*c^4 + 5*b^8*c^4 - 12*a^7*c^5 + 32*a^6*b*c^5 + 30*a^5*b^2*c^5 - 112*a^4*b^3*c^5 + 8*a^3*b^4*c^5 + 96*a^2*b^5*
c^5 - 26*a*b^6*c^5 - 16*b^7*c^5 - 104*a^6*c^6 + 192*a^5*b*c^6 + 149*a^4*b^2*c^6 - 336*a^3*b^3*c^6 - 30*a^2*b^4
*c^6 + 112*a*b^5*c^6 + 17*b^6*c^6 - 276*a^5*c^7 + 320*a^4*b*c^7 + 292*a^3*b^2*c^7 - 224*a^2*b^3*c^7 - 120*a*b^
4*c^7 - 304*a^4*c^8 + 128*a^3*b*c^8 + 237*a^2*b^2*c^8 + 24*a*b^3*c^8 - 17*b^4*c^8 - 116*a^3*c^9 - 96*a^2*b*c^9
+ 62*a*b^2*c^9 + 16*b^3*c^9 + 24*a^2*c^10 - 64*a*b*c^10 - 5*b^2*c^10 + 20*a*c^11)*abs(c)) - ((2*a^3*b^6 - 6*a
^2*b^7 + 6*a*b^8 - 2*b^9 - 18*a^4*b^4*c + 56*a^3*b^5*c - 54*a^2*b^6*c + 12*a*b^7*c + 4*b^8*c + 48*a^5*b^2*c^2
- 160*a^4*b^3*c^2 + 140*a^3*b^4*c^2 + 12*a^2*b^5*c^2 - 38*a*b^6*c^2 - 2*b^7*c^2 - 32*a^6*c^3 + 128*a^5*b*c^3 -
64*a^4*b^2*c^3 - 160*a^3*b^3*c^3 + 110*a^2*b^4*c^3 + 20*a*b^5*c^3 - 64*a^5*c^4 + 192*a^4*b*c^4 - 80*a^3*b^2*c
^4 - 64*a^2*b^3*c^4 - 32*a^4*c^5 + 64*a^3*b*c^5 + 3*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c)
)*sqrt(b^2 - 4*a*c)*a^3*b^4 - 2*(b^2 - 4*a*c)*a^3*b^4 - 5*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a -
b + c))*sqrt(b^2 - 4*a*c)*a^2*b^5 + 6*(b^2 - 4*a*c)*a^2*b^5 - 3*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)
*(a - b + c))*sqrt(b^2 - 4*a*c)*a*b^6 - 6*(b^2 - 4*a*c)*a*b^6 + 5*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*
c)*(a - b + c))*sqrt(b^2 - 4*a*c)*b^7 + 2*(b^2 - 4*a*c)*b^7 - 15*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c
)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^4*b^2*c + 10*(b^2 - 4*a*c)*a^4*b^2*c + 28*sqrt(a^2 - a*b + b*c - c^2 - sqrt
(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^3*b^3*c - 32*(b^2 - 4*a*c)*a^3*b^3*c + 27*sqrt(a^2 - a*b + b*c
- c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^2*b^4*c + 30*(b^2 - 4*a*c)*a^2*b^4*c - 38*sqrt(a^2
- a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a*b^5*c - 4*(b^2 - 4*a*c)*a*b^5*c - 6*sqr
t(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*b^6*c - 4*(b^2 - 4*a*c)*b^6*c + 12*
sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^5*c^2 - 8*(b^2 - 4*a*c)*a^5*c^
2 - 32*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^4*b*c^2 + 32*(b^2 - 4*a
*c)*a^4*b*c^2 - 74*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^3*b^2*c^2 -
20*(b^2 - 4*a*c)*a^3*b^2*c^2 + 94*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*
c)*a^2*b^3*c^2 - 28*(b^2 - 4*a*c)*a^2*b^3*c^2 + 31*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))
*sqrt(b^2 - 4*a*c)*a*b^4*c^2 + 22*(b^2 - 4*a*c)*a*b^4*c^2 + 5*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(
a - b + c))*sqrt(b^2 - 4*a*c)*b^5*c^2 + 2*(b^2 - 4*a*c)*b^5*c^2 + 56*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4
*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^4*c^3 - 16*(b^2 - 4*a*c)*a^4*c^3 - 88*sqrt(a^2 - a*b + b*c - c^2 - sqrt
(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^3*b*c^3 + 48*(b^2 - 4*a*c)*a^3*b*c^3 - 23*sqrt(a^2 - a*b + b*c
- c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^2*b^2*c^3 - 22*(b^2 - 4*a*c)*a^2*b^2*c^3 - 30*sqrt(
a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a*b^3*c^3 - 12*(b^2 - 4*a*c)*a*b^3*c^
3 - 20*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^3*c^4 - 8*(b^2 - 4*a*c)
*a^3*c^4 + 40*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^2*b*c^4 + 16*(b^
2 - 4*a*c)*a^2*b*c^4)*c^2*abs(a - b + c) - (4*a^3*b^6*c - 4*a^2*b^7*c - 4*a*b^8*c + 4*b^9*c - 36*a^4*b^4*c^2 +
32*a^3*b^5*c^2 + 52*a^2*b^6*c^2 - 40*a*b^7*c^2 - 8*b^8*c^2 + 96*a^5*b^2*c^3 - 64*a^4*b^3*c^3 - 232*a^3*b^4*c^
3 + 120*a^2*b^5*c^3 + 84*a*b^6*c^3 + 4*b^7*c^3 - 64*a^6*c^4 + 384*a^4*b^2*c^4 - 64*a^3*b^3*c^4 - 292*a^2*b^4*c
^4 - 40*a*b^5*c^4 - 128*a^5*c^5 - 128*a^4*b*c^5 + 352*a^3*b^2*c^5 + 128*a^2*b^3*c^5 - 64*a^4*c^6 - 128*a^3*b*c
^6 - 3*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*a^4*b^4*c + 2*sqrt(a^2 - a*b + b*c - c^2 -
sqrt(b^2 - 4*a*c)*(a - b + c))*a^3*b^5*c + 8*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*a^2*b
^6*c - 2*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*a*b^7*c - 5*sqrt(a^2 - a*b + b*c - c^2 -
sqrt(b^2 - 4*a*c)*(a - b + c))*b^8*c + 15*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*a^5*b^2*
c^2 - 7*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*a^4*b^3*c^2 - 62*sqrt(a^2 - a*b + b*c - c^
2 - sqrt(b^2 - 4*a*c)*(a - b + c))*a^3*b^4*c^2 + 51*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c)
)*a*b^6*c^2 + 11*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*b^7*c^2 - 12*sqrt(a^2 - a*b + b*c
- c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*a^6*c^3 - 4*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c)
)*a^5*b*c^3 + 137*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*a^4*b^2*c^3 + 56*sqrt(a^2 - a*b
+ b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*a^3*b^3*c^3 - 164*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*
(a - b + c))*a^2*b^4*c^3 - 86*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*a*b^5*c^3 - 11*sqrt(
a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*b^6*c^3 - 68*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*
a*c)*(a - b + c))*a^5*c^4 - 96*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*a^4*b*c^4 + 169*sqr
t(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*a^3*b^2*c^4 + 197*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b
^2 - 4*a*c)*(a - b + c))*a^2*b^3*c^4 + 71*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*a*b^4*c^
4 + 5*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*b^5*c^4 - 36*sqrt(a^2 - a*b + b*c - c^2 - sq
rt(b^2 - 4*a*c)*(a - b + c))*a^4*c^5 - 116*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*a^3*b*c
^5 - 113*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*a^2*b^2*c^5 - 30*sqrt(a^2 - a*b + b*c - c
^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*a*b^3*c^5 + 20*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))
*a^3*c^6 + 40*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*a^2*b*c^6 - 4*(b^2 - 4*a*c)*a^3*b^4*
c + 4*(b^2 - 4*a*c)*a^2*b^5*c + 4*(b^2 - 4*a*c)*a*b^6*c - 4*(b^2 - 4*a*c)*b^7*c + 20*(b^2 - 4*a*c)*a^4*b^2*c^2
- 16*(b^2 - 4*a*c)*a^3*b^3*c^2 - 36*(b^2 - 4*a*c)*a^2*b^4*c^2 + 24*(b^2 - 4*a*c)*a*b^5*c^2 + 8*(b^2 - 4*a*c)*
b^6*c^2 - 16*(b^2 - 4*a*c)*a^5*c^3 + 88*(b^2 - 4*a*c)*a^3*b^2*c^3 - 24*(b^2 - 4*a*c)*a^2*b^3*c^3 - 52*(b^2 - 4
*a*c)*a*b^4*c^3 - 4*(b^2 - 4*a*c)*b^5*c^3 - 32*(b^2 - 4*a*c)*a^4*c^4 - 32*(b^2 - 4*a*c)*a^3*b*c^4 + 84*(b^2 -
4*a*c)*a^2*b^2*c^4 + 24*(b^2 - 4*a*c)*a*b^3*c^4 - 16*(b^2 - 4*a*c)*a^3*c^5 - 32*(b^2 - 4*a*c)*a^2*b*c^5)*abs(a
- b + c)*abs(c) + (2*a^4*b^5*c^2 - 6*a^3*b^6*c^2 + 6*a^2*b^7*c^2 - 2*a*b^8*c^2 - 14*a^5*b^3*c^3 + 44*a^4*b^4*
c^3 - 44*a^3*b^5*c^3 + 14*a^2*b^6*c^3 - 2*a*b^7*c^3 + 2*b^8*c^3 + 24*a^6*b*c^4 - 84*a^5*b^2*c^4 + 82*a^4*b^3*c
^4 - 20*a^3*b^4*c^4 + 12*a^2*b^5*c^4 - 10*a*b^6*c^4 - 4*b^7*c^4 + 16*a^6*c^5 - 8*a^5*b*c^5 - 20*a^4*b^2*c^5 -
10*a^3*b^3*c^5 - 8*a^2*b^4*c^5 + 30*a*b^5*c^5 + 2*b^6*c^5 + 16*a^5*c^6 - 24*a^4*b*c^6 + 68*a^3*b^2*c^6 - 58*a^
2*b^3*c^6 - 16*a*b^4*c^6 - 16*a^4*c^7 + 8*a^3*b*c^7 + 36*a^2*b^2*c^7 - 16*a^3*c^8 + 3*sqrt(a^2 - a*b + b*c - c
^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^4*b^3*c^2 - 2*(b^2 - 4*a*c)*a^4*b^3*c^2 - 5*sqrt(a^2 -
a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^3*b^4*c^2 + 6*(b^2 - 4*a*c)*a^3*b^4*c^2
- 3*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^2*b^5*c^2 - 6*(b^2 - 4*a*c
)*a^2*b^5*c^2 + 5*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a*b^6*c^2 + 2*
(b^2 - 4*a*c)*a*b^6*c^2 - 9*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^5*
b*c^3 + 6*(b^2 - 4*a*c)*a^5*b*c^3 + 18*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 -
4*a*c)*a^4*b^2*c^3 - 20*(b^2 - 4*a*c)*a^4*b^2*c^3 + 18*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b +
c))*sqrt(b^2 - 4*a*c)*a^3*b^3*c^3 + 20*(b^2 - 4*a*c)*a^3*b^3*c^3 - 23*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 -
4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^2*b^4*c^3 - 6*(b^2 - 4*a*c)*a^2*b^4*c^3 - 3*sqrt(a^2 - a*b + b*c - c^
2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a*b^5*c^3 + 2*(b^2 - 4*a*c)*a*b^5*c^3 - 5*sqrt(a^2 - a*b
+ b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*b^6*c^3 - 2*(b^2 - 4*a*c)*b^6*c^3 - 6*sqrt(a^2
- a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^5*c^4 + 4*(b^2 - 4*a*c)*a^5*c^4 - 29*sq
rt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^4*b*c^4 - 2*(b^2 - 4*a*c)*a^4*b*
c^4 + 30*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^3*b^2*c^4 - 4*(b^2 -
4*a*c)*a^3*b^2*c^4 - 2*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^2*b^3*c
^4 - 4*(b^2 - 4*a*c)*a^2*b^3*c^4 + 33*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4
*a*c)*a*b^4*c^4 + 2*(b^2 - 4*a*c)*a*b^4*c^4 + 6*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sq
rt(b^2 - 4*a*c)*b^5*c^4 + 4*(b^2 - 4*a*c)*b^5*c^4 - 22*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b +
c))*sqrt(b^2 - 4*a*c)*a^4*c^5 + 4*(b^2 - 4*a*c)*a^4*c^5 + 41*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(
a - b + c))*sqrt(b^2 - 4*a*c)*a^3*b*c^5 - 6*(b^2 - 4*a*c)*a^3*b*c^5 - 68*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2
- 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^2*b^2*c^5 + 16*(b^2 - 4*a*c)*a^2*b^2*c^5 - 19*sqrt(a^2 - a*b + b*c
- c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a*b^3*c^5 - 14*(b^2 - 4*a*c)*a*b^3*c^5 - 5*sqrt(a^2 -
a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*b^4*c^5 - 2*(b^2 - 4*a*c)*b^4*c^5 + 38*sqr
t(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^3*c^6 - 4*(b^2 - 4*a*c)*a^3*c^6 -
3*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^2*b*c^6 + 2*(b^2 - 4*a*c)*a
^2*b*c^6 + 20*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a*b^2*c^6 + 8*(b^2
- 4*a*c)*a*b^2*c^6 - 10*sqrt(a^2 - a*b + b*c - c^2 - sqrt(b^2 - 4*a*c)*(a - b + c))*sqrt(b^2 - 4*a*c)*a^2*c^7
- 4*(b^2 - 4*a*c)*a^2*c^7)*abs(a - b + c))*(pi*floor(1/2*x/pi + 1/2) + arctan(2*sqrt(1/2)*tan(1/2*x)/sqrt((2*
a*c^3 - 2*c^4 - sqrt(-4*(a*c^3 + b*c^3 + c^4)*(a*c^3 - b*c^3 + c^4) + 4*(a*c^3 - c^4)^2))/(a*c^3 - b*c^3 + c^4
))))/((3*a^6*b^2*c^4 - 8*a^5*b^3*c^4 - a^4*b^4*c^4 + 16*a^3*b^5*c^4 - 7*a^2*b^6*c^4 - 8*a*b^7*c^4 + 5*b^8*c^4
- 12*a^7*c^5 + 32*a^6*b*c^5 + 30*a^5*b^2*c^5 - 112*a^4*b^3*c^5 + 8*a^3*b^4*c^5 + 96*a^2*b^5*c^5 - 26*a*b^6*c^5
- 16*b^7*c^5 - 104*a^6*c^6 + 192*a^5*b*c^6 + 149*a^4*b^2*c^6 - 336*a^3*b^3*c^6 - 30*a^2*b^4*c^6 + 112*a*b^5*c
^6 + 17*b^6*c^6 - 276*a^5*c^7 + 320*a^4*b*c^7 + 292*a^3*b^2*c^7 - 224*a^2*b^3*c^7 - 120*a*b^4*c^7 - 304*a^4*c^
8 + 128*a^3*b*c^8 + 237*a^2*b^2*c^8 + 24*a*b^3*c^8 - 17*b^4*c^8 - 116*a^3*c^9 - 96*a^2*b*c^9 + 62*a*b^2*c^9 +
16*b^3*c^9 + 24*a^2*c^10 - 64*a*b*c^10 - 5*b^2*c^10 + 20*a*c^11)*abs(c)) + 1/2*(2*b^2 - 2*a*c + c^2)*x/c^3 - (
2*b*tan(1/2*x)^3 + c*tan(1/2*x)^3 + 2*b*tan(1/2*x) - c*tan(1/2*x))/((tan(1/2*x)^2 + 1)^2*c^2)

Mupad [B] (verification not implemented)

Time = 16.24 (sec) , antiderivative size = 45364, normalized size of antiderivative = 139.15 \[ \int \frac {\cos ^4(x)}{a+b \cos (x)+c \cos ^2(x)} \, dx=\text {Too large to display} \]

[In]

int(cos(x)^4/(a + b*cos(x) + c*cos(x)^2),x)

[Out]

atan(((((2048*(12*a^3*c^11 - 28*a^4*c^10 - 44*a^5*c^9 + 72*a^6*c^8 + 88*a^7*c^7 + 12*a^8*c^6 + b^5*c^9 - 4*b^6
*c^8 + 10*b^7*c^7 - 20*b^8*c^6 + 29*b^9*c^5 - 30*b^10*c^4 + 26*b^11*c^3 - 12*b^12*c^2 - 6*a*b^3*c^10 + 27*a*b^
4*c^9 - 72*a*b^5*c^8 + 154*a*b^6*c^7 - 238*a*b^7*c^6 + 251*a*b^8*c^5 - 228*a*b^9*c^4 + 98*a*b^10*c^3 + 20*a*b^
11*c^2 + 8*a^2*b*c^11 - 68*a^3*b*c^10 + 112*a^4*b*c^9 + 100*a^5*b*c^8 - 200*a^6*b*c^7 - 96*a^7*b*c^6 - 47*a^2*
b^2*c^10 + 145*a^2*b^3*c^9 - 354*a^2*b^4*c^8 + 612*a^2*b^5*c^7 - 655*a^2*b^6*c^6 + 635*a^2*b^7*c^5 - 202*a^2*b
^8*c^4 - 222*a^2*b^9*c^3 + 4*a^2*b^10*c^2 + 239*a^3*b^2*c^9 - 524*a^3*b^3*c^8 + 536*a^3*b^4*c^7 - 564*a^3*b^5*
c^6 - 115*a^3*b^6*c^5 + 856*a^3*b^7*c^4 + 2*a^3*b^8*c^3 - 20*a^3*b^9*c^2 - 37*a^4*b^2*c^8 + 9*a^4*b^3*c^7 + 58
3*a^4*b^4*c^6 - 1362*a^4*b^5*c^5 - 152*a^4*b^6*c^4 + 156*a^4*b^7*c^3 + 8*a^4*b^8*c^2 - 399*a^5*b^2*c^7 + 904*a
^5*b^3*c^6 + 394*a^5*b^4*c^5 - 388*a^5*b^5*c^4 - 60*a^5*b^6*c^3 - 340*a^6*b^2*c^6 + 364*a^6*b^3*c^5 + 136*a^6*
b^4*c^4 - 100*a^7*b^2*c^5))/c^8 + (((2048*(16*a^3*c^13 - 32*a^2*c^14 + 176*a^4*c^12 + 176*a^5*c^11 + 48*a^6*c^
10 - 2*b^4*c^12 + 6*b^5*c^11 - 18*b^6*c^10 + 26*b^7*c^9 - 12*b^8*c^8 + 16*a*b^2*c^13 - 40*a*b^3*c^12 + 122*a*b
^4*c^11 - 192*a*b^5*c^10 + 74*a*b^6*c^9 + 20*a*b^7*c^8 + 64*a^2*b*c^13 - 144*a^3*b*c^12 - 352*a^4*b*c^11 - 144
*a^5*b*c^10 - 204*a^2*b^2*c^12 + 388*a^2*b^3*c^11 - 50*a^2*b^4*c^10 - 182*a^2*b^5*c^9 + 4*a^2*b^6*c^8 - 260*a^
3*b^2*c^11 + 496*a^3*b^3*c^10 + 10*a^3*b^4*c^9 - 20*a^3*b^5*c^8 - 148*a^4*b^2*c^10 + 116*a^4*b^3*c^9 + 8*a^4*b
^4*c^8 - 44*a^5*b^2*c^9))/c^8 - (2048*tan(x/2)*(-(a^2*b^8 - b^10 + 8*a^5*c^5 + 8*a^6*c^4 - b^7*(-(4*a*c - b^2)
^3)^(1/2) - 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) - 52*a^2*b^6*c^2 + 96*a^3*b^4*c^3 - 66*a^4*b^2*c^4
+ 33*a^4*b^4*c^2 - 38*a^5*b^2*c^3 + 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c
- b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(
-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^
4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2)*(32*a*c^16 - 64*a^2*c^15 - 128*a^3*c^14 + 64*a^4
*c^13 + 96*a^5*c^12 - 8*b^2*c^15 + 24*b^3*c^14 - 32*b^4*c^13 + 32*b^5*c^12 - 24*b^6*c^11 + 8*b^7*c^10 + 144*a*
b^2*c^14 - 200*a*b^3*c^13 + 184*a*b^4*c^12 - 56*a*b^5*c^11 - 8*a*b^6*c^10 + 288*a^2*b*c^14 + 352*a^3*b*c^13 -
32*a^4*b*c^12 - 320*a^2*b^2*c^13 + 8*a^2*b^3*c^12 + 96*a^2*b^4*c^11 - 8*a^2*b^5*c^10 - 272*a^3*b^2*c^12 + 40*a
^3*b^3*c^11 + 8*a^3*b^4*c^10 - 56*a^4*b^2*c^11 - 96*a*b*c^15))/c^8)*(-(a^2*b^8 - b^10 + 8*a^5*c^5 + 8*a^6*c^4
- b^7*(-(4*a*c - b^2)^3)^(1/2) - 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) - 52*a^2*b^6*c^2 + 96*a^3*b^4
*c^3 - 66*a^4*b^2*c^4 + 33*a^4*b^4*c^2 - 38*a^5*b^2*c^3 + 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) -
4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3
)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 -
8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2) - (2048*tan(x/2)*(8*a*c^14
- 64*a^2*c^13 + 80*a^3*c^12 + 168*a^4*c^11 - 192*a^5*c^10 - 136*a^6*c^9 + 72*a^7*c^8 - 2*b^2*c^13 + 6*b^3*c^12
- 17*b^4*c^11 + 33*b^5*c^10 - 49*b^6*c^9 + 61*b^7*c^8 - 52*b^8*c^7 + 36*b^9*c^6 - 24*b^10*c^5 + 8*b^11*c^4 +
84*a*b^2*c^12 - 178*a*b^3*c^11 + 295*a*b^4*c^10 - 416*a*b^5*c^9 + 375*a*b^6*c^8 - 308*a*b^7*c^7 + 244*a*b^8*c^
6 - 72*a*b^9*c^5 - 8*a*b^10*c^4 + 184*a^2*b*c^12 - 328*a^3*b*c^11 - 16*a^4*b*c^10 + 496*a^5*b*c^9 - 88*a^6*b*c
^8 - 416*a^2*b^2*c^11 + 770*a^2*b^3*c^10 - 723*a^2*b^4*c^9 + 779*a^2*b^5*c^8 - 732*a^2*b^6*c^7 + 80*a^2*b^7*c^
6 + 112*a^2*b^8*c^5 - 8*a^2*b^9*c^4 + 180*a^3*b^2*c^10 - 494*a^3*b^3*c^9 + 521*a^3*b^4*c^8 + 572*a^3*b^5*c^7 -
424*a^3*b^6*c^6 + 56*a^3*b^7*c^5 + 8*a^3*b^8*c^4 + 234*a^4*b^2*c^9 - 1152*a^4*b^3*c^8 + 416*a^4*b^4*c^7 - 140
*a^4*b^5*c^6 - 72*a^4*b^6*c^5 + 64*a^5*b^2*c^8 + 192*a^5*b^3*c^7 + 220*a^5*b^4*c^6 - 256*a^6*b^2*c^7 - 24*a*b*
c^13))/c^8)*(-(a^2*b^8 - b^10 + 8*a^5*c^5 + 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) - 10*a^3*b^6*c + a^2*b^5*
(-(4*a*c - b^2)^3)^(1/2) - 52*a^2*b^6*c^2 + 96*a^3*b^4*c^3 - 66*a^4*b^2*c^4 + 33*a^4*b^4*c^2 - 38*a^5*b^2*c^3
+ 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*
a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2
*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^
6 - 8*a^3*b^2*c^7)))^(1/2))*(-(a^2*b^8 - b^10 + 8*a^5*c^5 + 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) - 10*a^3*
b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) - 52*a^2*b^6*c^2 + 96*a^3*b^4*c^3 - 66*a^4*b^2*c^4 + 33*a^4*b^4*c^2 -
38*a^5*b^2*c^3 + 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3
*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(
1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*
c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2) - (2048*tan(x/2)*(20*a*b^12 + 4*b^12*c - 4*b^13 - 40*a^2*b^11 + 40*
a^3*b^10 - 20*a^4*b^9 + 4*a^5*b^8 + 2*a^4*c^9 - 18*a^5*c^8 + 38*a^6*c^7 + 2*a^7*c^6 - 44*a^8*c^5 + 12*a^9*c^4
+ b^8*c^5 - b^9*c^4 + 4*b^10*c^3 - 4*b^11*c^2 - 8*a*b^6*c^6 + 4*a*b^7*c^5 - 31*a*b^8*c^4 + 20*a*b^9*c^3 - 20*a
*b^10*c^2 - 160*a^2*b^10*c + 320*a^3*b^9*c + 26*a^4*b*c^8 - 300*a^4*b^8*c - 84*a^5*b*c^7 + 136*a^5*b^7*c + 2*a
^6*b*c^6 - 24*a^6*b^6*c + 168*a^7*b*c^5 - 92*a^8*b*c^4 + 20*a^2*b^4*c^7 + 8*a^2*b^5*c^6 + 82*a^2*b^6*c^5 + 6*a
^2*b^7*c^4 + 8*a^2*b^8*c^3 - 44*a^2*b^9*c^2 - 16*a^3*b^2*c^8 - 40*a^3*b^3*c^7 - 104*a^3*b^4*c^6 - 132*a^3*b^5*
c^5 + 34*a^3*b^6*c^4 + 72*a^3*b^7*c^3 + 460*a^3*b^8*c^2 + 82*a^4*b^2*c^7 + 174*a^4*b^3*c^6 + 41*a^4*b^4*c^5 -
149*a^4*b^5*c^4 - 660*a^4*b^6*c^3 - 900*a^4*b^7*c^2 - 90*a^5*b^2*c^6 + 96*a^5*b^3*c^5 + 541*a^5*b^4*c^4 + 1156
*a^5*b^5*c^3 + 764*a^5*b^6*c^2 - 204*a^6*b^2*c^5 - 704*a^6*b^3*c^4 - 840*a^6*b^4*c^3 - 300*a^6*b^5*c^2 + 384*a
^7*b^2*c^4 + 272*a^7*b^3*c^3 + 44*a^7*b^4*c^2 - 32*a^8*b^2*c^3 + 24*a*b^11*c))/c^8)*(-(a^2*b^8 - b^10 + 8*a^5*
c^5 + 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) - 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) - 52*a^2*b^6*
c^2 + 96*a^3*b^4*c^3 - 66*a^4*b^2*c^4 + 33*a^4*b^4*c^2 - 38*a^5*b^2*c^3 + 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c -
b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(
-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*
c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2)*1i - (((204
8*(12*a^3*c^11 - 28*a^4*c^10 - 44*a^5*c^9 + 72*a^6*c^8 + 88*a^7*c^7 + 12*a^8*c^6 + b^5*c^9 - 4*b^6*c^8 + 10*b^
7*c^7 - 20*b^8*c^6 + 29*b^9*c^5 - 30*b^10*c^4 + 26*b^11*c^3 - 12*b^12*c^2 - 6*a*b^3*c^10 + 27*a*b^4*c^9 - 72*a
*b^5*c^8 + 154*a*b^6*c^7 - 238*a*b^7*c^6 + 251*a*b^8*c^5 - 228*a*b^9*c^4 + 98*a*b^10*c^3 + 20*a*b^11*c^2 + 8*a
^2*b*c^11 - 68*a^3*b*c^10 + 112*a^4*b*c^9 + 100*a^5*b*c^8 - 200*a^6*b*c^7 - 96*a^7*b*c^6 - 47*a^2*b^2*c^10 + 1
45*a^2*b^3*c^9 - 354*a^2*b^4*c^8 + 612*a^2*b^5*c^7 - 655*a^2*b^6*c^6 + 635*a^2*b^7*c^5 - 202*a^2*b^8*c^4 - 222
*a^2*b^9*c^3 + 4*a^2*b^10*c^2 + 239*a^3*b^2*c^9 - 524*a^3*b^3*c^8 + 536*a^3*b^4*c^7 - 564*a^3*b^5*c^6 - 115*a^
3*b^6*c^5 + 856*a^3*b^7*c^4 + 2*a^3*b^8*c^3 - 20*a^3*b^9*c^2 - 37*a^4*b^2*c^8 + 9*a^4*b^3*c^7 + 583*a^4*b^4*c^
6 - 1362*a^4*b^5*c^5 - 152*a^4*b^6*c^4 + 156*a^4*b^7*c^3 + 8*a^4*b^8*c^2 - 399*a^5*b^2*c^7 + 904*a^5*b^3*c^6 +
394*a^5*b^4*c^5 - 388*a^5*b^5*c^4 - 60*a^5*b^6*c^3 - 340*a^6*b^2*c^6 + 364*a^6*b^3*c^5 + 136*a^6*b^4*c^4 - 10
0*a^7*b^2*c^5))/c^8 + (((2048*(16*a^3*c^13 - 32*a^2*c^14 + 176*a^4*c^12 + 176*a^5*c^11 + 48*a^6*c^10 - 2*b^4*c
^12 + 6*b^5*c^11 - 18*b^6*c^10 + 26*b^7*c^9 - 12*b^8*c^8 + 16*a*b^2*c^13 - 40*a*b^3*c^12 + 122*a*b^4*c^11 - 19
2*a*b^5*c^10 + 74*a*b^6*c^9 + 20*a*b^7*c^8 + 64*a^2*b*c^13 - 144*a^3*b*c^12 - 352*a^4*b*c^11 - 144*a^5*b*c^10
- 204*a^2*b^2*c^12 + 388*a^2*b^3*c^11 - 50*a^2*b^4*c^10 - 182*a^2*b^5*c^9 + 4*a^2*b^6*c^8 - 260*a^3*b^2*c^11 +
496*a^3*b^3*c^10 + 10*a^3*b^4*c^9 - 20*a^3*b^5*c^8 - 148*a^4*b^2*c^10 + 116*a^4*b^3*c^9 + 8*a^4*b^4*c^8 - 44*
a^5*b^2*c^9))/c^8 + (2048*tan(x/2)*(-(a^2*b^8 - b^10 + 8*a^5*c^5 + 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) -
10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) - 52*a^2*b^6*c^2 + 96*a^3*b^4*c^3 - 66*a^4*b^2*c^4 + 33*a^4*b^
4*c^2 - 38*a^5*b^2*c^3 + 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1
/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^
2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a
^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2)*(32*a*c^16 - 64*a^2*c^15 - 128*a^3*c^14 + 64*a^4*c^13 + 96*a
^5*c^12 - 8*b^2*c^15 + 24*b^3*c^14 - 32*b^4*c^13 + 32*b^5*c^12 - 24*b^6*c^11 + 8*b^7*c^10 + 144*a*b^2*c^14 - 2
00*a*b^3*c^13 + 184*a*b^4*c^12 - 56*a*b^5*c^11 - 8*a*b^6*c^10 + 288*a^2*b*c^14 + 352*a^3*b*c^13 - 32*a^4*b*c^1
2 - 320*a^2*b^2*c^13 + 8*a^2*b^3*c^12 + 96*a^2*b^4*c^11 - 8*a^2*b^5*c^10 - 272*a^3*b^2*c^12 + 40*a^3*b^3*c^11
+ 8*a^3*b^4*c^10 - 56*a^4*b^2*c^11 - 96*a*b*c^15))/c^8)*(-(a^2*b^8 - b^10 + 8*a^5*c^5 + 8*a^6*c^4 - b^7*(-(4*a
*c - b^2)^3)^(1/2) - 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) - 52*a^2*b^6*c^2 + 96*a^3*b^4*c^3 - 66*a^
4*b^2*c^4 + 33*a^4*b^4*c^2 - 38*a^5*b^2*c^3 + 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*
(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*
a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9
+ 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2) + (2048*tan(x/2)*(8*a*c^14 - 64*a^2*c^1
3 + 80*a^3*c^12 + 168*a^4*c^11 - 192*a^5*c^10 - 136*a^6*c^9 + 72*a^7*c^8 - 2*b^2*c^13 + 6*b^3*c^12 - 17*b^4*c^
11 + 33*b^5*c^10 - 49*b^6*c^9 + 61*b^7*c^8 - 52*b^8*c^7 + 36*b^9*c^6 - 24*b^10*c^5 + 8*b^11*c^4 + 84*a*b^2*c^1
2 - 178*a*b^3*c^11 + 295*a*b^4*c^10 - 416*a*b^5*c^9 + 375*a*b^6*c^8 - 308*a*b^7*c^7 + 244*a*b^8*c^6 - 72*a*b^9
*c^5 - 8*a*b^10*c^4 + 184*a^2*b*c^12 - 328*a^3*b*c^11 - 16*a^4*b*c^10 + 496*a^5*b*c^9 - 88*a^6*b*c^8 - 416*a^2
*b^2*c^11 + 770*a^2*b^3*c^10 - 723*a^2*b^4*c^9 + 779*a^2*b^5*c^8 - 732*a^2*b^6*c^7 + 80*a^2*b^7*c^6 + 112*a^2*
b^8*c^5 - 8*a^2*b^9*c^4 + 180*a^3*b^2*c^10 - 494*a^3*b^3*c^9 + 521*a^3*b^4*c^8 + 572*a^3*b^5*c^7 - 424*a^3*b^6
*c^6 + 56*a^3*b^7*c^5 + 8*a^3*b^8*c^4 + 234*a^4*b^2*c^9 - 1152*a^4*b^3*c^8 + 416*a^4*b^4*c^7 - 140*a^4*b^5*c^6
- 72*a^4*b^6*c^5 + 64*a^5*b^2*c^8 + 192*a^5*b^3*c^7 + 220*a^5*b^4*c^6 - 256*a^6*b^2*c^7 - 24*a*b*c^13))/c^8)*
(-(a^2*b^8 - b^10 + 8*a^5*c^5 + 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) - 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b
^2)^3)^(1/2) - 52*a^2*b^6*c^2 + 96*a^3*b^4*c^3 - 66*a^4*b^2*c^4 + 33*a^4*b^4*c^2 - 38*a^5*b^2*c^3 + 12*a*b^8*c
+ 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3
)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a
^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^
2*c^7)))^(1/2))*(-(a^2*b^8 - b^10 + 8*a^5*c^5 + 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) - 10*a^3*b^6*c + a^2*
b^5*(-(4*a*c - b^2)^3)^(1/2) - 52*a^2*b^6*c^2 + 96*a^3*b^4*c^3 - 66*a^4*b^2*c^4 + 33*a^4*b^4*c^2 - 38*a^5*b^2*
c^3 + 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(
-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16
*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^
4*c^6 - 8*a^3*b^2*c^7)))^(1/2) + (2048*tan(x/2)*(20*a*b^12 + 4*b^12*c - 4*b^13 - 40*a^2*b^11 + 40*a^3*b^10 - 2
0*a^4*b^9 + 4*a^5*b^8 + 2*a^4*c^9 - 18*a^5*c^8 + 38*a^6*c^7 + 2*a^7*c^6 - 44*a^8*c^5 + 12*a^9*c^4 + b^8*c^5 -
b^9*c^4 + 4*b^10*c^3 - 4*b^11*c^2 - 8*a*b^6*c^6 + 4*a*b^7*c^5 - 31*a*b^8*c^4 + 20*a*b^9*c^3 - 20*a*b^10*c^2 -
160*a^2*b^10*c + 320*a^3*b^9*c + 26*a^4*b*c^8 - 300*a^4*b^8*c - 84*a^5*b*c^7 + 136*a^5*b^7*c + 2*a^6*b*c^6 - 2
4*a^6*b^6*c + 168*a^7*b*c^5 - 92*a^8*b*c^4 + 20*a^2*b^4*c^7 + 8*a^2*b^5*c^6 + 82*a^2*b^6*c^5 + 6*a^2*b^7*c^4 +
8*a^2*b^8*c^3 - 44*a^2*b^9*c^2 - 16*a^3*b^2*c^8 - 40*a^3*b^3*c^7 - 104*a^3*b^4*c^6 - 132*a^3*b^5*c^5 + 34*a^3
*b^6*c^4 + 72*a^3*b^7*c^3 + 460*a^3*b^8*c^2 + 82*a^4*b^2*c^7 + 174*a^4*b^3*c^6 + 41*a^4*b^4*c^5 - 149*a^4*b^5*
c^4 - 660*a^4*b^6*c^3 - 900*a^4*b^7*c^2 - 90*a^5*b^2*c^6 + 96*a^5*b^3*c^5 + 541*a^5*b^4*c^4 + 1156*a^5*b^5*c^3
+ 764*a^5*b^6*c^2 - 204*a^6*b^2*c^5 - 704*a^6*b^3*c^4 - 840*a^6*b^4*c^3 - 300*a^6*b^5*c^2 + 384*a^7*b^2*c^4 +
272*a^7*b^3*c^3 + 44*a^7*b^4*c^2 - 32*a^8*b^2*c^3 + 24*a*b^11*c))/c^8)*(-(a^2*b^8 - b^10 + 8*a^5*c^5 + 8*a^6*
c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) - 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) - 52*a^2*b^6*c^2 + 96*a^3
*b^4*c^3 - 66*a^4*b^2*c^4 + 33*a^4*b^4*c^2 - 38*a^5*b^2*c^3 + 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2
) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^
2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^
6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2)*1i)/((4096*(16*a^5*b^7
- 4*a^4*b^8 - 24*a^6*b^6 + 16*a^7*b^5 - 4*a^8*b^4 + 3*a^6*c^6 - 10*a^7*c^5 + a^8*c^4 + 14*a^9*c^3 + 4*a^4*b^7*
c - 2*a^5*b*c^6 + 4*a^5*b^6*c + 6*a^6*b*c^5 - 40*a^6*b^5*c + 4*a^7*b*c^4 + 56*a^7*b^4*c - 22*a^8*b*c^3 - 28*a^
8*b^3*c + 12*a^9*b*c^2 + 4*a^9*b^2*c + a^4*b^3*c^5 - a^4*b^4*c^4 + 4*a^4*b^5*c^3 - 4*a^4*b^6*c^2 - a^5*b^2*c^5
- 8*a^5*b^3*c^4 + 10*a^6*b^2*c^4 - 4*a^6*b^3*c^3 - 8*a^6*b^4*c^2 + 4*a^7*b^2*c^3 + 48*a^7*b^3*c^2 - 48*a^8*b^
2*c^2))/c^8 + (((2048*(12*a^3*c^11 - 28*a^4*c^10 - 44*a^5*c^9 + 72*a^6*c^8 + 88*a^7*c^7 + 12*a^8*c^6 + b^5*c^9
- 4*b^6*c^8 + 10*b^7*c^7 - 20*b^8*c^6 + 29*b^9*c^5 - 30*b^10*c^4 + 26*b^11*c^3 - 12*b^12*c^2 - 6*a*b^3*c^10 +
27*a*b^4*c^9 - 72*a*b^5*c^8 + 154*a*b^6*c^7 - 238*a*b^7*c^6 + 251*a*b^8*c^5 - 228*a*b^9*c^4 + 98*a*b^10*c^3 +
20*a*b^11*c^2 + 8*a^2*b*c^11 - 68*a^3*b*c^10 + 112*a^4*b*c^9 + 100*a^5*b*c^8 - 200*a^6*b*c^7 - 96*a^7*b*c^6 -
47*a^2*b^2*c^10 + 145*a^2*b^3*c^9 - 354*a^2*b^4*c^8 + 612*a^2*b^5*c^7 - 655*a^2*b^6*c^6 + 635*a^2*b^7*c^5 - 2
02*a^2*b^8*c^4 - 222*a^2*b^9*c^3 + 4*a^2*b^10*c^2 + 239*a^3*b^2*c^9 - 524*a^3*b^3*c^8 + 536*a^3*b^4*c^7 - 564*
a^3*b^5*c^6 - 115*a^3*b^6*c^5 + 856*a^3*b^7*c^4 + 2*a^3*b^8*c^3 - 20*a^3*b^9*c^2 - 37*a^4*b^2*c^8 + 9*a^4*b^3*
c^7 + 583*a^4*b^4*c^6 - 1362*a^4*b^5*c^5 - 152*a^4*b^6*c^4 + 156*a^4*b^7*c^3 + 8*a^4*b^8*c^2 - 399*a^5*b^2*c^7
+ 904*a^5*b^3*c^6 + 394*a^5*b^4*c^5 - 388*a^5*b^5*c^4 - 60*a^5*b^6*c^3 - 340*a^6*b^2*c^6 + 364*a^6*b^3*c^5 +
136*a^6*b^4*c^4 - 100*a^7*b^2*c^5))/c^8 + (((2048*(16*a^3*c^13 - 32*a^2*c^14 + 176*a^4*c^12 + 176*a^5*c^11 + 4
8*a^6*c^10 - 2*b^4*c^12 + 6*b^5*c^11 - 18*b^6*c^10 + 26*b^7*c^9 - 12*b^8*c^8 + 16*a*b^2*c^13 - 40*a*b^3*c^12 +
122*a*b^4*c^11 - 192*a*b^5*c^10 + 74*a*b^6*c^9 + 20*a*b^7*c^8 + 64*a^2*b*c^13 - 144*a^3*b*c^12 - 352*a^4*b*c^
11 - 144*a^5*b*c^10 - 204*a^2*b^2*c^12 + 388*a^2*b^3*c^11 - 50*a^2*b^4*c^10 - 182*a^2*b^5*c^9 + 4*a^2*b^6*c^8
- 260*a^3*b^2*c^11 + 496*a^3*b^3*c^10 + 10*a^3*b^4*c^9 - 20*a^3*b^5*c^8 - 148*a^4*b^2*c^10 + 116*a^4*b^3*c^9 +
8*a^4*b^4*c^8 - 44*a^5*b^2*c^9))/c^8 - (2048*tan(x/2)*(-(a^2*b^8 - b^10 + 8*a^5*c^5 + 8*a^6*c^4 - b^7*(-(4*a*
c - b^2)^3)^(1/2) - 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) - 52*a^2*b^6*c^2 + 96*a^3*b^4*c^3 - 66*a^4
*b^2*c^4 + 33*a^4*b^4*c^2 - 38*a^5*b^2*c^3 + 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(
-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a
*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 +
10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2)*(32*a*c^16 - 64*a^2*c^15 - 128*a^3*c^14
+ 64*a^4*c^13 + 96*a^5*c^12 - 8*b^2*c^15 + 24*b^3*c^14 - 32*b^4*c^13 + 32*b^5*c^12 - 24*b^6*c^11 + 8*b^7*c^10
+ 144*a*b^2*c^14 - 200*a*b^3*c^13 + 184*a*b^4*c^12 - 56*a*b^5*c^11 - 8*a*b^6*c^10 + 288*a^2*b*c^14 + 352*a^3*b
*c^13 - 32*a^4*b*c^12 - 320*a^2*b^2*c^13 + 8*a^2*b^3*c^12 + 96*a^2*b^4*c^11 - 8*a^2*b^5*c^10 - 272*a^3*b^2*c^1
2 + 40*a^3*b^3*c^11 + 8*a^3*b^4*c^10 - 56*a^4*b^2*c^11 - 96*a*b*c^15))/c^8)*(-(a^2*b^8 - b^10 + 8*a^5*c^5 + 8*
a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) - 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) - 52*a^2*b^6*c^2 + 96
*a^3*b^4*c^3 - 66*a^4*b^2*c^4 + 33*a^4*b^4*c^2 - 38*a^5*b^2*c^3 + 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^
(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c
- b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^
6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2) - (2048*tan(x/2)*(8
*a*c^14 - 64*a^2*c^13 + 80*a^3*c^12 + 168*a^4*c^11 - 192*a^5*c^10 - 136*a^6*c^9 + 72*a^7*c^8 - 2*b^2*c^13 + 6*
b^3*c^12 - 17*b^4*c^11 + 33*b^5*c^10 - 49*b^6*c^9 + 61*b^7*c^8 - 52*b^8*c^7 + 36*b^9*c^6 - 24*b^10*c^5 + 8*b^1
1*c^4 + 84*a*b^2*c^12 - 178*a*b^3*c^11 + 295*a*b^4*c^10 - 416*a*b^5*c^9 + 375*a*b^6*c^8 - 308*a*b^7*c^7 + 244*
a*b^8*c^6 - 72*a*b^9*c^5 - 8*a*b^10*c^4 + 184*a^2*b*c^12 - 328*a^3*b*c^11 - 16*a^4*b*c^10 + 496*a^5*b*c^9 - 88
*a^6*b*c^8 - 416*a^2*b^2*c^11 + 770*a^2*b^3*c^10 - 723*a^2*b^4*c^9 + 779*a^2*b^5*c^8 - 732*a^2*b^6*c^7 + 80*a^
2*b^7*c^6 + 112*a^2*b^8*c^5 - 8*a^2*b^9*c^4 + 180*a^3*b^2*c^10 - 494*a^3*b^3*c^9 + 521*a^3*b^4*c^8 + 572*a^3*b
^5*c^7 - 424*a^3*b^6*c^6 + 56*a^3*b^7*c^5 + 8*a^3*b^8*c^4 + 234*a^4*b^2*c^9 - 1152*a^4*b^3*c^8 + 416*a^4*b^4*c
^7 - 140*a^4*b^5*c^6 - 72*a^4*b^6*c^5 + 64*a^5*b^2*c^8 + 192*a^5*b^3*c^7 + 220*a^5*b^4*c^6 - 256*a^6*b^2*c^7 -
24*a*b*c^13))/c^8)*(-(a^2*b^8 - b^10 + 8*a^5*c^5 + 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) - 10*a^3*b^6*c +
a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) - 52*a^2*b^6*c^2 + 96*a^3*b^4*c^3 - 66*a^4*b^2*c^4 + 33*a^4*b^4*c^2 - 38*a^5*
b^2*c^3 + 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c
^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2
*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^
2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2))*(-(a^2*b^8 - b^10 + 8*a^5*c^5 + 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) -
10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) - 52*a^2*b^6*c^2 + 96*a^3*b^4*c^3 - 66*a^4*b^2*c^4 + 33*a^4*b
^4*c^2 - 38*a^5*b^2*c^3 + 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(
1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b
^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*
a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2) - (2048*tan(x/2)*(20*a*b^12 + 4*b^12*c - 4*b^13 - 40*a^2*b^
11 + 40*a^3*b^10 - 20*a^4*b^9 + 4*a^5*b^8 + 2*a^4*c^9 - 18*a^5*c^8 + 38*a^6*c^7 + 2*a^7*c^6 - 44*a^8*c^5 + 12*
a^9*c^4 + b^8*c^5 - b^9*c^4 + 4*b^10*c^3 - 4*b^11*c^2 - 8*a*b^6*c^6 + 4*a*b^7*c^5 - 31*a*b^8*c^4 + 20*a*b^9*c^
3 - 20*a*b^10*c^2 - 160*a^2*b^10*c + 320*a^3*b^9*c + 26*a^4*b*c^8 - 300*a^4*b^8*c - 84*a^5*b*c^7 + 136*a^5*b^7
*c + 2*a^6*b*c^6 - 24*a^6*b^6*c + 168*a^7*b*c^5 - 92*a^8*b*c^4 + 20*a^2*b^4*c^7 + 8*a^2*b^5*c^6 + 82*a^2*b^6*c
^5 + 6*a^2*b^7*c^4 + 8*a^2*b^8*c^3 - 44*a^2*b^9*c^2 - 16*a^3*b^2*c^8 - 40*a^3*b^3*c^7 - 104*a^3*b^4*c^6 - 132*
a^3*b^5*c^5 + 34*a^3*b^6*c^4 + 72*a^3*b^7*c^3 + 460*a^3*b^8*c^2 + 82*a^4*b^2*c^7 + 174*a^4*b^3*c^6 + 41*a^4*b^
4*c^5 - 149*a^4*b^5*c^4 - 660*a^4*b^6*c^3 - 900*a^4*b^7*c^2 - 90*a^5*b^2*c^6 + 96*a^5*b^3*c^5 + 541*a^5*b^4*c^
4 + 1156*a^5*b^5*c^3 + 764*a^5*b^6*c^2 - 204*a^6*b^2*c^5 - 704*a^6*b^3*c^4 - 840*a^6*b^4*c^3 - 300*a^6*b^5*c^2
+ 384*a^7*b^2*c^4 + 272*a^7*b^3*c^3 + 44*a^7*b^4*c^2 - 32*a^8*b^2*c^3 + 24*a*b^11*c))/c^8)*(-(a^2*b^8 - b^10
+ 8*a^5*c^5 + 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) - 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) - 52*
a^2*b^6*c^2 + 96*a^3*b^4*c^3 - 66*a^4*b^2*c^4 + 33*a^4*b^4*c^2 - 38*a^5*b^2*c^3 + 12*a*b^8*c + 4*a^3*b*c^3*(-(
4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b
^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^
8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2) + (
((2048*(12*a^3*c^11 - 28*a^4*c^10 - 44*a^5*c^9 + 72*a^6*c^8 + 88*a^7*c^7 + 12*a^8*c^6 + b^5*c^9 - 4*b^6*c^8 +
10*b^7*c^7 - 20*b^8*c^6 + 29*b^9*c^5 - 30*b^10*c^4 + 26*b^11*c^3 - 12*b^12*c^2 - 6*a*b^3*c^10 + 27*a*b^4*c^9 -
72*a*b^5*c^8 + 154*a*b^6*c^7 - 238*a*b^7*c^6 + 251*a*b^8*c^5 - 228*a*b^9*c^4 + 98*a*b^10*c^3 + 20*a*b^11*c^2
+ 8*a^2*b*c^11 - 68*a^3*b*c^10 + 112*a^4*b*c^9 + 100*a^5*b*c^8 - 200*a^6*b*c^7 - 96*a^7*b*c^6 - 47*a^2*b^2*c^1
0 + 145*a^2*b^3*c^9 - 354*a^2*b^4*c^8 + 612*a^2*b^5*c^7 - 655*a^2*b^6*c^6 + 635*a^2*b^7*c^5 - 202*a^2*b^8*c^4
- 222*a^2*b^9*c^3 + 4*a^2*b^10*c^2 + 239*a^3*b^2*c^9 - 524*a^3*b^3*c^8 + 536*a^3*b^4*c^7 - 564*a^3*b^5*c^6 - 1
15*a^3*b^6*c^5 + 856*a^3*b^7*c^4 + 2*a^3*b^8*c^3 - 20*a^3*b^9*c^2 - 37*a^4*b^2*c^8 + 9*a^4*b^3*c^7 + 583*a^4*b
^4*c^6 - 1362*a^4*b^5*c^5 - 152*a^4*b^6*c^4 + 156*a^4*b^7*c^3 + 8*a^4*b^8*c^2 - 399*a^5*b^2*c^7 + 904*a^5*b^3*
c^6 + 394*a^5*b^4*c^5 - 388*a^5*b^5*c^4 - 60*a^5*b^6*c^3 - 340*a^6*b^2*c^6 + 364*a^6*b^3*c^5 + 136*a^6*b^4*c^4
- 100*a^7*b^2*c^5))/c^8 + (((2048*(16*a^3*c^13 - 32*a^2*c^14 + 176*a^4*c^12 + 176*a^5*c^11 + 48*a^6*c^10 - 2*
b^4*c^12 + 6*b^5*c^11 - 18*b^6*c^10 + 26*b^7*c^9 - 12*b^8*c^8 + 16*a*b^2*c^13 - 40*a*b^3*c^12 + 122*a*b^4*c^11
- 192*a*b^5*c^10 + 74*a*b^6*c^9 + 20*a*b^7*c^8 + 64*a^2*b*c^13 - 144*a^3*b*c^12 - 352*a^4*b*c^11 - 144*a^5*b*
c^10 - 204*a^2*b^2*c^12 + 388*a^2*b^3*c^11 - 50*a^2*b^4*c^10 - 182*a^2*b^5*c^9 + 4*a^2*b^6*c^8 - 260*a^3*b^2*c
^11 + 496*a^3*b^3*c^10 + 10*a^3*b^4*c^9 - 20*a^3*b^5*c^8 - 148*a^4*b^2*c^10 + 116*a^4*b^3*c^9 + 8*a^4*b^4*c^8
- 44*a^5*b^2*c^9))/c^8 + (2048*tan(x/2)*(-(a^2*b^8 - b^10 + 8*a^5*c^5 + 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/
2) - 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) - 52*a^2*b^6*c^2 + 96*a^3*b^4*c^3 - 66*a^4*b^2*c^4 + 33*a
^4*b^4*c^2 - 38*a^5*b^2*c^3 + 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^
3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c
- b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 -
32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2)*(32*a*c^16 - 64*a^2*c^15 - 128*a^3*c^14 + 64*a^4*c^13 +
96*a^5*c^12 - 8*b^2*c^15 + 24*b^3*c^14 - 32*b^4*c^13 + 32*b^5*c^12 - 24*b^6*c^11 + 8*b^7*c^10 + 144*a*b^2*c^1
4 - 200*a*b^3*c^13 + 184*a*b^4*c^12 - 56*a*b^5*c^11 - 8*a*b^6*c^10 + 288*a^2*b*c^14 + 352*a^3*b*c^13 - 32*a^4*
b*c^12 - 320*a^2*b^2*c^13 + 8*a^2*b^3*c^12 + 96*a^2*b^4*c^11 - 8*a^2*b^5*c^10 - 272*a^3*b^2*c^12 + 40*a^3*b^3*
c^11 + 8*a^3*b^4*c^10 - 56*a^4*b^2*c^11 - 96*a*b*c^15))/c^8)*(-(a^2*b^8 - b^10 + 8*a^5*c^5 + 8*a^6*c^4 - b^7*(
-(4*a*c - b^2)^3)^(1/2) - 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) - 52*a^2*b^6*c^2 + 96*a^3*b^4*c^3 -
66*a^4*b^2*c^4 + 33*a^4*b^4*c^2 - 38*a^5*b^2*c^3 + 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b
^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2)
+ 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2
*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2) + (2048*tan(x/2)*(8*a*c^14 - 64*a^
2*c^13 + 80*a^3*c^12 + 168*a^4*c^11 - 192*a^5*c^10 - 136*a^6*c^9 + 72*a^7*c^8 - 2*b^2*c^13 + 6*b^3*c^12 - 17*b
^4*c^11 + 33*b^5*c^10 - 49*b^6*c^9 + 61*b^7*c^8 - 52*b^8*c^7 + 36*b^9*c^6 - 24*b^10*c^5 + 8*b^11*c^4 + 84*a*b^
2*c^12 - 178*a*b^3*c^11 + 295*a*b^4*c^10 - 416*a*b^5*c^9 + 375*a*b^6*c^8 - 308*a*b^7*c^7 + 244*a*b^8*c^6 - 72*
a*b^9*c^5 - 8*a*b^10*c^4 + 184*a^2*b*c^12 - 328*a^3*b*c^11 - 16*a^4*b*c^10 + 496*a^5*b*c^9 - 88*a^6*b*c^8 - 41
6*a^2*b^2*c^11 + 770*a^2*b^3*c^10 - 723*a^2*b^4*c^9 + 779*a^2*b^5*c^8 - 732*a^2*b^6*c^7 + 80*a^2*b^7*c^6 + 112
*a^2*b^8*c^5 - 8*a^2*b^9*c^4 + 180*a^3*b^2*c^10 - 494*a^3*b^3*c^9 + 521*a^3*b^4*c^8 + 572*a^3*b^5*c^7 - 424*a^
3*b^6*c^6 + 56*a^3*b^7*c^5 + 8*a^3*b^8*c^4 + 234*a^4*b^2*c^9 - 1152*a^4*b^3*c^8 + 416*a^4*b^4*c^7 - 140*a^4*b^
5*c^6 - 72*a^4*b^6*c^5 + 64*a^5*b^2*c^8 + 192*a^5*b^3*c^7 + 220*a^5*b^4*c^6 - 256*a^6*b^2*c^7 - 24*a*b*c^13))/
c^8)*(-(a^2*b^8 - b^10 + 8*a^5*c^5 + 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) - 10*a^3*b^6*c + a^2*b^5*(-(4*a*
c - b^2)^3)^(1/2) - 52*a^2*b^6*c^2 + 96*a^3*b^4*c^3 - 66*a^4*b^2*c^4 + 33*a^4*b^4*c^2 - 38*a^5*b^2*c^3 + 12*a*
b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b
^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 +
32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a
^3*b^2*c^7)))^(1/2))*(-(a^2*b^8 - b^10 + 8*a^5*c^5 + 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) - 10*a^3*b^6*c +
a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) - 52*a^2*b^6*c^2 + 96*a^3*b^4*c^3 - 66*a^4*b^2*c^4 + 33*a^4*b^4*c^2 - 38*a^5
*b^2*c^3 + 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*
c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(
2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a
^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2) + (2048*tan(x/2)*(20*a*b^12 + 4*b^12*c - 4*b^13 - 40*a^2*b^11 + 40*a^3*b^1
0 - 20*a^4*b^9 + 4*a^5*b^8 + 2*a^4*c^9 - 18*a^5*c^8 + 38*a^6*c^7 + 2*a^7*c^6 - 44*a^8*c^5 + 12*a^9*c^4 + b^8*c
^5 - b^9*c^4 + 4*b^10*c^3 - 4*b^11*c^2 - 8*a*b^6*c^6 + 4*a*b^7*c^5 - 31*a*b^8*c^4 + 20*a*b^9*c^3 - 20*a*b^10*c
^2 - 160*a^2*b^10*c + 320*a^3*b^9*c + 26*a^4*b*c^8 - 300*a^4*b^8*c - 84*a^5*b*c^7 + 136*a^5*b^7*c + 2*a^6*b*c^
6 - 24*a^6*b^6*c + 168*a^7*b*c^5 - 92*a^8*b*c^4 + 20*a^2*b^4*c^7 + 8*a^2*b^5*c^6 + 82*a^2*b^6*c^5 + 6*a^2*b^7*
c^4 + 8*a^2*b^8*c^3 - 44*a^2*b^9*c^2 - 16*a^3*b^2*c^8 - 40*a^3*b^3*c^7 - 104*a^3*b^4*c^6 - 132*a^3*b^5*c^5 + 3
4*a^3*b^6*c^4 + 72*a^3*b^7*c^3 + 460*a^3*b^8*c^2 + 82*a^4*b^2*c^7 + 174*a^4*b^3*c^6 + 41*a^4*b^4*c^5 - 149*a^4
*b^5*c^4 - 660*a^4*b^6*c^3 - 900*a^4*b^7*c^2 - 90*a^5*b^2*c^6 + 96*a^5*b^3*c^5 + 541*a^5*b^4*c^4 + 1156*a^5*b^
5*c^3 + 764*a^5*b^6*c^2 - 204*a^6*b^2*c^5 - 704*a^6*b^3*c^4 - 840*a^6*b^4*c^3 - 300*a^6*b^5*c^2 + 384*a^7*b^2*
c^4 + 272*a^7*b^3*c^3 + 44*a^7*b^4*c^2 - 32*a^8*b^2*c^3 + 24*a*b^11*c))/c^8)*(-(a^2*b^8 - b^10 + 8*a^5*c^5 + 8
*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) - 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) - 52*a^2*b^6*c^2 + 9
6*a^3*b^4*c^3 - 66*a^4*b^2*c^4 + 33*a^4*b^4*c^2 - 38*a^5*b^2*c^3 + 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)
^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c
- b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b
^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2)))*(-(a^2*b^8 - b^1
0 + 8*a^5*c^5 + 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) - 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) - 5
2*a^2*b^6*c^2 + 96*a^3*b^4*c^3 - 66*a^4*b^2*c^4 + 33*a^4*b^4*c^2 - 38*a^5*b^2*c^3 + 12*a*b^8*c + 4*a^3*b*c^3*(
-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2
*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*
c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2)*2
i - ((tan(x/2)*(2*b - c))/c^2 + (tan(x/2)^3*(2*b + c))/c^2)/(2*tan(x/2)^2 + tan(x/2)^4 + 1) + atan(((((2048*(1
2*a^3*c^11 - 28*a^4*c^10 - 44*a^5*c^9 + 72*a^6*c^8 + 88*a^7*c^7 + 12*a^8*c^6 + b^5*c^9 - 4*b^6*c^8 + 10*b^7*c^
7 - 20*b^8*c^6 + 29*b^9*c^5 - 30*b^10*c^4 + 26*b^11*c^3 - 12*b^12*c^2 - 6*a*b^3*c^10 + 27*a*b^4*c^9 - 72*a*b^5
*c^8 + 154*a*b^6*c^7 - 238*a*b^7*c^6 + 251*a*b^8*c^5 - 228*a*b^9*c^4 + 98*a*b^10*c^3 + 20*a*b^11*c^2 + 8*a^2*b
*c^11 - 68*a^3*b*c^10 + 112*a^4*b*c^9 + 100*a^5*b*c^8 - 200*a^6*b*c^7 - 96*a^7*b*c^6 - 47*a^2*b^2*c^10 + 145*a
^2*b^3*c^9 - 354*a^2*b^4*c^8 + 612*a^2*b^5*c^7 - 655*a^2*b^6*c^6 + 635*a^2*b^7*c^5 - 202*a^2*b^8*c^4 - 222*a^2
*b^9*c^3 + 4*a^2*b^10*c^2 + 239*a^3*b^2*c^9 - 524*a^3*b^3*c^8 + 536*a^3*b^4*c^7 - 564*a^3*b^5*c^6 - 115*a^3*b^
6*c^5 + 856*a^3*b^7*c^4 + 2*a^3*b^8*c^3 - 20*a^3*b^9*c^2 - 37*a^4*b^2*c^8 + 9*a^4*b^3*c^7 + 583*a^4*b^4*c^6 -
1362*a^4*b^5*c^5 - 152*a^4*b^6*c^4 + 156*a^4*b^7*c^3 + 8*a^4*b^8*c^2 - 399*a^5*b^2*c^7 + 904*a^5*b^3*c^6 + 394
*a^5*b^4*c^5 - 388*a^5*b^5*c^4 - 60*a^5*b^6*c^3 - 340*a^6*b^2*c^6 + 364*a^6*b^3*c^5 + 136*a^6*b^4*c^4 - 100*a^
7*b^2*c^5))/c^8 + (((2048*(16*a^3*c^13 - 32*a^2*c^14 + 176*a^4*c^12 + 176*a^5*c^11 + 48*a^6*c^10 - 2*b^4*c^12
+ 6*b^5*c^11 - 18*b^6*c^10 + 26*b^7*c^9 - 12*b^8*c^8 + 16*a*b^2*c^13 - 40*a*b^3*c^12 + 122*a*b^4*c^11 - 192*a*
b^5*c^10 + 74*a*b^6*c^9 + 20*a*b^7*c^8 + 64*a^2*b*c^13 - 144*a^3*b*c^12 - 352*a^4*b*c^11 - 144*a^5*b*c^10 - 20
4*a^2*b^2*c^12 + 388*a^2*b^3*c^11 - 50*a^2*b^4*c^10 - 182*a^2*b^5*c^9 + 4*a^2*b^6*c^8 - 260*a^3*b^2*c^11 + 496
*a^3*b^3*c^10 + 10*a^3*b^4*c^9 - 20*a^3*b^5*c^8 - 148*a^4*b^2*c^10 + 116*a^4*b^3*c^9 + 8*a^4*b^4*c^8 - 44*a^5*
b^2*c^9))/c^8 - (2048*tan(x/2)*((b^10 - a^2*b^8 - 8*a^5*c^5 - 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) + 10*a^
3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) + 52*a^2*b^6*c^2 - 96*a^3*b^4*c^3 + 66*a^4*b^2*c^4 - 33*a^4*b^4*c^2
+ 38*a^5*b^2*c^3 - 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) +
3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)
^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^
2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2)*(32*a*c^16 - 64*a^2*c^15 - 128*a^3*c^14 + 64*a^4*c^13 + 96*a^5*c^
12 - 8*b^2*c^15 + 24*b^3*c^14 - 32*b^4*c^13 + 32*b^5*c^12 - 24*b^6*c^11 + 8*b^7*c^10 + 144*a*b^2*c^14 - 200*a*
b^3*c^13 + 184*a*b^4*c^12 - 56*a*b^5*c^11 - 8*a*b^6*c^10 + 288*a^2*b*c^14 + 352*a^3*b*c^13 - 32*a^4*b*c^12 - 3
20*a^2*b^2*c^13 + 8*a^2*b^3*c^12 + 96*a^2*b^4*c^11 - 8*a^2*b^5*c^10 - 272*a^3*b^2*c^12 + 40*a^3*b^3*c^11 + 8*a
^3*b^4*c^10 - 56*a^4*b^2*c^11 - 96*a*b*c^15))/c^8)*((b^10 - a^2*b^8 - 8*a^5*c^5 - 8*a^6*c^4 - b^7*(-(4*a*c - b
^2)^3)^(1/2) + 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) + 52*a^2*b^6*c^2 - 96*a^3*b^4*c^3 + 66*a^4*b^2*
c^4 - 33*a^4*b^4*c^2 + 38*a^5*b^2*c^3 - 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a
*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*
c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a
*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2) - (2048*tan(x/2)*(8*a*c^14 - 64*a^2*c^13 + 80
*a^3*c^12 + 168*a^4*c^11 - 192*a^5*c^10 - 136*a^6*c^9 + 72*a^7*c^8 - 2*b^2*c^13 + 6*b^3*c^12 - 17*b^4*c^11 + 3
3*b^5*c^10 - 49*b^6*c^9 + 61*b^7*c^8 - 52*b^8*c^7 + 36*b^9*c^6 - 24*b^10*c^5 + 8*b^11*c^4 + 84*a*b^2*c^12 - 17
8*a*b^3*c^11 + 295*a*b^4*c^10 - 416*a*b^5*c^9 + 375*a*b^6*c^8 - 308*a*b^7*c^7 + 244*a*b^8*c^6 - 72*a*b^9*c^5 -
8*a*b^10*c^4 + 184*a^2*b*c^12 - 328*a^3*b*c^11 - 16*a^4*b*c^10 + 496*a^5*b*c^9 - 88*a^6*b*c^8 - 416*a^2*b^2*c
^11 + 770*a^2*b^3*c^10 - 723*a^2*b^4*c^9 + 779*a^2*b^5*c^8 - 732*a^2*b^6*c^7 + 80*a^2*b^7*c^6 + 112*a^2*b^8*c^
5 - 8*a^2*b^9*c^4 + 180*a^3*b^2*c^10 - 494*a^3*b^3*c^9 + 521*a^3*b^4*c^8 + 572*a^3*b^5*c^7 - 424*a^3*b^6*c^6 +
56*a^3*b^7*c^5 + 8*a^3*b^8*c^4 + 234*a^4*b^2*c^9 - 1152*a^4*b^3*c^8 + 416*a^4*b^4*c^7 - 140*a^4*b^5*c^6 - 72*
a^4*b^6*c^5 + 64*a^5*b^2*c^8 + 192*a^5*b^3*c^7 + 220*a^5*b^4*c^6 - 256*a^6*b^2*c^7 - 24*a*b*c^13))/c^8)*((b^10
- a^2*b^8 - 8*a^5*c^5 - 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) + 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^
(1/2) + 52*a^2*b^6*c^2 - 96*a^3*b^4*c^3 + 66*a^4*b^2*c^4 - 33*a^4*b^4*c^2 + 38*a^5*b^2*c^3 - 12*a*b^8*c + 4*a^
3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2)
- 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9
+ 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7))
)^(1/2))*((b^10 - a^2*b^8 - 8*a^5*c^5 - 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) + 10*a^3*b^6*c + a^2*b^5*(-(4
*a*c - b^2)^3)^(1/2) + 52*a^2*b^6*c^2 - 96*a^3*b^4*c^3 + 66*a^4*b^2*c^4 - 33*a^4*b^4*c^2 + 38*a^5*b^2*c^3 - 12
*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c
- b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^1
0 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 -
8*a^3*b^2*c^7)))^(1/2) - (2048*tan(x/2)*(20*a*b^12 + 4*b^12*c - 4*b^13 - 40*a^2*b^11 + 40*a^3*b^10 - 20*a^4*b^
9 + 4*a^5*b^8 + 2*a^4*c^9 - 18*a^5*c^8 + 38*a^6*c^7 + 2*a^7*c^6 - 44*a^8*c^5 + 12*a^9*c^4 + b^8*c^5 - b^9*c^4
+ 4*b^10*c^3 - 4*b^11*c^2 - 8*a*b^6*c^6 + 4*a*b^7*c^5 - 31*a*b^8*c^4 + 20*a*b^9*c^3 - 20*a*b^10*c^2 - 160*a^2*
b^10*c + 320*a^3*b^9*c + 26*a^4*b*c^8 - 300*a^4*b^8*c - 84*a^5*b*c^7 + 136*a^5*b^7*c + 2*a^6*b*c^6 - 24*a^6*b^
6*c + 168*a^7*b*c^5 - 92*a^8*b*c^4 + 20*a^2*b^4*c^7 + 8*a^2*b^5*c^6 + 82*a^2*b^6*c^5 + 6*a^2*b^7*c^4 + 8*a^2*b
^8*c^3 - 44*a^2*b^9*c^2 - 16*a^3*b^2*c^8 - 40*a^3*b^3*c^7 - 104*a^3*b^4*c^6 - 132*a^3*b^5*c^5 + 34*a^3*b^6*c^4
+ 72*a^3*b^7*c^3 + 460*a^3*b^8*c^2 + 82*a^4*b^2*c^7 + 174*a^4*b^3*c^6 + 41*a^4*b^4*c^5 - 149*a^4*b^5*c^4 - 66
0*a^4*b^6*c^3 - 900*a^4*b^7*c^2 - 90*a^5*b^2*c^6 + 96*a^5*b^3*c^5 + 541*a^5*b^4*c^4 + 1156*a^5*b^5*c^3 + 764*a
^5*b^6*c^2 - 204*a^6*b^2*c^5 - 704*a^6*b^3*c^4 - 840*a^6*b^4*c^3 - 300*a^6*b^5*c^2 + 384*a^7*b^2*c^4 + 272*a^7
*b^3*c^3 + 44*a^7*b^4*c^2 - 32*a^8*b^2*c^3 + 24*a*b^11*c))/c^8)*((b^10 - a^2*b^8 - 8*a^5*c^5 - 8*a^6*c^4 - b^7
*(-(4*a*c - b^2)^3)^(1/2) + 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) + 52*a^2*b^6*c^2 - 96*a^3*b^4*c^3
+ 66*a^4*b^2*c^4 - 33*a^4*b^4*c^2 + 38*a^5*b^2*c^3 - 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3
*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/
2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b
^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2)*1i - (((2048*(12*a^3*c^11 - 28*a
^4*c^10 - 44*a^5*c^9 + 72*a^6*c^8 + 88*a^7*c^7 + 12*a^8*c^6 + b^5*c^9 - 4*b^6*c^8 + 10*b^7*c^7 - 20*b^8*c^6 +
29*b^9*c^5 - 30*b^10*c^4 + 26*b^11*c^3 - 12*b^12*c^2 - 6*a*b^3*c^10 + 27*a*b^4*c^9 - 72*a*b^5*c^8 + 154*a*b^6*
c^7 - 238*a*b^7*c^6 + 251*a*b^8*c^5 - 228*a*b^9*c^4 + 98*a*b^10*c^3 + 20*a*b^11*c^2 + 8*a^2*b*c^11 - 68*a^3*b*
c^10 + 112*a^4*b*c^9 + 100*a^5*b*c^8 - 200*a^6*b*c^7 - 96*a^7*b*c^6 - 47*a^2*b^2*c^10 + 145*a^2*b^3*c^9 - 354*
a^2*b^4*c^8 + 612*a^2*b^5*c^7 - 655*a^2*b^6*c^6 + 635*a^2*b^7*c^5 - 202*a^2*b^8*c^4 - 222*a^2*b^9*c^3 + 4*a^2*
b^10*c^2 + 239*a^3*b^2*c^9 - 524*a^3*b^3*c^8 + 536*a^3*b^4*c^7 - 564*a^3*b^5*c^6 - 115*a^3*b^6*c^5 + 856*a^3*b
^7*c^4 + 2*a^3*b^8*c^3 - 20*a^3*b^9*c^2 - 37*a^4*b^2*c^8 + 9*a^4*b^3*c^7 + 583*a^4*b^4*c^6 - 1362*a^4*b^5*c^5
- 152*a^4*b^6*c^4 + 156*a^4*b^7*c^3 + 8*a^4*b^8*c^2 - 399*a^5*b^2*c^7 + 904*a^5*b^3*c^6 + 394*a^5*b^4*c^5 - 38
8*a^5*b^5*c^4 - 60*a^5*b^6*c^3 - 340*a^6*b^2*c^6 + 364*a^6*b^3*c^5 + 136*a^6*b^4*c^4 - 100*a^7*b^2*c^5))/c^8 +
(((2048*(16*a^3*c^13 - 32*a^2*c^14 + 176*a^4*c^12 + 176*a^5*c^11 + 48*a^6*c^10 - 2*b^4*c^12 + 6*b^5*c^11 - 18
*b^6*c^10 + 26*b^7*c^9 - 12*b^8*c^8 + 16*a*b^2*c^13 - 40*a*b^3*c^12 + 122*a*b^4*c^11 - 192*a*b^5*c^10 + 74*a*b
^6*c^9 + 20*a*b^7*c^8 + 64*a^2*b*c^13 - 144*a^3*b*c^12 - 352*a^4*b*c^11 - 144*a^5*b*c^10 - 204*a^2*b^2*c^12 +
388*a^2*b^3*c^11 - 50*a^2*b^4*c^10 - 182*a^2*b^5*c^9 + 4*a^2*b^6*c^8 - 260*a^3*b^2*c^11 + 496*a^3*b^3*c^10 + 1
0*a^3*b^4*c^9 - 20*a^3*b^5*c^8 - 148*a^4*b^2*c^10 + 116*a^4*b^3*c^9 + 8*a^4*b^4*c^8 - 44*a^5*b^2*c^9))/c^8 + (
2048*tan(x/2)*((b^10 - a^2*b^8 - 8*a^5*c^5 - 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) + 10*a^3*b^6*c + a^2*b^5
*(-(4*a*c - b^2)^3)^(1/2) + 52*a^2*b^6*c^2 - 96*a^3*b^4*c^3 + 66*a^4*b^2*c^4 - 33*a^4*b^4*c^2 + 38*a^5*b^2*c^3
- 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4
*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^
2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c
^6 - 8*a^3*b^2*c^7)))^(1/2)*(32*a*c^16 - 64*a^2*c^15 - 128*a^3*c^14 + 64*a^4*c^13 + 96*a^5*c^12 - 8*b^2*c^15 +
24*b^3*c^14 - 32*b^4*c^13 + 32*b^5*c^12 - 24*b^6*c^11 + 8*b^7*c^10 + 144*a*b^2*c^14 - 200*a*b^3*c^13 + 184*a*
b^4*c^12 - 56*a*b^5*c^11 - 8*a*b^6*c^10 + 288*a^2*b*c^14 + 352*a^3*b*c^13 - 32*a^4*b*c^12 - 320*a^2*b^2*c^13 +
8*a^2*b^3*c^12 + 96*a^2*b^4*c^11 - 8*a^2*b^5*c^10 - 272*a^3*b^2*c^12 + 40*a^3*b^3*c^11 + 8*a^3*b^4*c^10 - 56*
a^4*b^2*c^11 - 96*a*b*c^15))/c^8)*((b^10 - a^2*b^8 - 8*a^5*c^5 - 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) + 10
*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) + 52*a^2*b^6*c^2 - 96*a^3*b^4*c^3 + 66*a^4*b^2*c^4 - 33*a^4*b^4*
c^2 + 38*a^5*b^2*c^3 - 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2
) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)
^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2
*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2) + (2048*tan(x/2)*(8*a*c^14 - 64*a^2*c^13 + 80*a^3*c^12 + 168*a
^4*c^11 - 192*a^5*c^10 - 136*a^6*c^9 + 72*a^7*c^8 - 2*b^2*c^13 + 6*b^3*c^12 - 17*b^4*c^11 + 33*b^5*c^10 - 49*b
^6*c^9 + 61*b^7*c^8 - 52*b^8*c^7 + 36*b^9*c^6 - 24*b^10*c^5 + 8*b^11*c^4 + 84*a*b^2*c^12 - 178*a*b^3*c^11 + 29
5*a*b^4*c^10 - 416*a*b^5*c^9 + 375*a*b^6*c^8 - 308*a*b^7*c^7 + 244*a*b^8*c^6 - 72*a*b^9*c^5 - 8*a*b^10*c^4 + 1
84*a^2*b*c^12 - 328*a^3*b*c^11 - 16*a^4*b*c^10 + 496*a^5*b*c^9 - 88*a^6*b*c^8 - 416*a^2*b^2*c^11 + 770*a^2*b^3
*c^10 - 723*a^2*b^4*c^9 + 779*a^2*b^5*c^8 - 732*a^2*b^6*c^7 + 80*a^2*b^7*c^6 + 112*a^2*b^8*c^5 - 8*a^2*b^9*c^4
+ 180*a^3*b^2*c^10 - 494*a^3*b^3*c^9 + 521*a^3*b^4*c^8 + 572*a^3*b^5*c^7 - 424*a^3*b^6*c^6 + 56*a^3*b^7*c^5 +
8*a^3*b^8*c^4 + 234*a^4*b^2*c^9 - 1152*a^4*b^3*c^8 + 416*a^4*b^4*c^7 - 140*a^4*b^5*c^6 - 72*a^4*b^6*c^5 + 64*
a^5*b^2*c^8 + 192*a^5*b^3*c^7 + 220*a^5*b^4*c^6 - 256*a^6*b^2*c^7 - 24*a*b*c^13))/c^8)*((b^10 - a^2*b^8 - 8*a^
5*c^5 - 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) + 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) + 52*a^2*b^
6*c^2 - 96*a^3*b^4*c^3 + 66*a^4*b^2*c^4 - 33*a^4*b^4*c^2 + 38*a^5*b^2*c^3 - 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c
- b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2
*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^
4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2))*((b^10 -
a^2*b^8 - 8*a^5*c^5 - 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) + 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1
/2) + 52*a^2*b^6*c^2 - 96*a^3*b^4*c^3 + 66*a^4*b^2*c^4 - 33*a^4*b^4*c^2 + 38*a^5*b^2*c^3 - 12*a*b^8*c + 4*a^3*
b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) -
10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 +
16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^
(1/2) + (2048*tan(x/2)*(20*a*b^12 + 4*b^12*c - 4*b^13 - 40*a^2*b^11 + 40*a^3*b^10 - 20*a^4*b^9 + 4*a^5*b^8 + 2
*a^4*c^9 - 18*a^5*c^8 + 38*a^6*c^7 + 2*a^7*c^6 - 44*a^8*c^5 + 12*a^9*c^4 + b^8*c^5 - b^9*c^4 + 4*b^10*c^3 - 4*
b^11*c^2 - 8*a*b^6*c^6 + 4*a*b^7*c^5 - 31*a*b^8*c^4 + 20*a*b^9*c^3 - 20*a*b^10*c^2 - 160*a^2*b^10*c + 320*a^3*
b^9*c + 26*a^4*b*c^8 - 300*a^4*b^8*c - 84*a^5*b*c^7 + 136*a^5*b^7*c + 2*a^6*b*c^6 - 24*a^6*b^6*c + 168*a^7*b*c
^5 - 92*a^8*b*c^4 + 20*a^2*b^4*c^7 + 8*a^2*b^5*c^6 + 82*a^2*b^6*c^5 + 6*a^2*b^7*c^4 + 8*a^2*b^8*c^3 - 44*a^2*b
^9*c^2 - 16*a^3*b^2*c^8 - 40*a^3*b^3*c^7 - 104*a^3*b^4*c^6 - 132*a^3*b^5*c^5 + 34*a^3*b^6*c^4 + 72*a^3*b^7*c^3
+ 460*a^3*b^8*c^2 + 82*a^4*b^2*c^7 + 174*a^4*b^3*c^6 + 41*a^4*b^4*c^5 - 149*a^4*b^5*c^4 - 660*a^4*b^6*c^3 - 9
00*a^4*b^7*c^2 - 90*a^5*b^2*c^6 + 96*a^5*b^3*c^5 + 541*a^5*b^4*c^4 + 1156*a^5*b^5*c^3 + 764*a^5*b^6*c^2 - 204*
a^6*b^2*c^5 - 704*a^6*b^3*c^4 - 840*a^6*b^4*c^3 - 300*a^6*b^5*c^2 + 384*a^7*b^2*c^4 + 272*a^7*b^3*c^3 + 44*a^7
*b^4*c^2 - 32*a^8*b^2*c^3 + 24*a*b^11*c))/c^8)*((b^10 - a^2*b^8 - 8*a^5*c^5 - 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^
3)^(1/2) + 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) + 52*a^2*b^6*c^2 - 96*a^3*b^4*c^3 + 66*a^4*b^2*c^4
- 33*a^4*b^4*c^2 + 38*a^5*b^2*c^3 - 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c -
b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-
(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4
*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2)*1i)/((4096*(16*a^5*b^7 - 4*a^4*b^8 - 24*a^6*b^6 +
16*a^7*b^5 - 4*a^8*b^4 + 3*a^6*c^6 - 10*a^7*c^5 + a^8*c^4 + 14*a^9*c^3 + 4*a^4*b^7*c - 2*a^5*b*c^6 + 4*a^5*b^
6*c + 6*a^6*b*c^5 - 40*a^6*b^5*c + 4*a^7*b*c^4 + 56*a^7*b^4*c - 22*a^8*b*c^3 - 28*a^8*b^3*c + 12*a^9*b*c^2 + 4
*a^9*b^2*c + a^4*b^3*c^5 - a^4*b^4*c^4 + 4*a^4*b^5*c^3 - 4*a^4*b^6*c^2 - a^5*b^2*c^5 - 8*a^5*b^3*c^4 + 10*a^6*
b^2*c^4 - 4*a^6*b^3*c^3 - 8*a^6*b^4*c^2 + 4*a^7*b^2*c^3 + 48*a^7*b^3*c^2 - 48*a^8*b^2*c^2))/c^8 + (((2048*(12*
a^3*c^11 - 28*a^4*c^10 - 44*a^5*c^9 + 72*a^6*c^8 + 88*a^7*c^7 + 12*a^8*c^6 + b^5*c^9 - 4*b^6*c^8 + 10*b^7*c^7
- 20*b^8*c^6 + 29*b^9*c^5 - 30*b^10*c^4 + 26*b^11*c^3 - 12*b^12*c^2 - 6*a*b^3*c^10 + 27*a*b^4*c^9 - 72*a*b^5*c
^8 + 154*a*b^6*c^7 - 238*a*b^7*c^6 + 251*a*b^8*c^5 - 228*a*b^9*c^4 + 98*a*b^10*c^3 + 20*a*b^11*c^2 + 8*a^2*b*c
^11 - 68*a^3*b*c^10 + 112*a^4*b*c^9 + 100*a^5*b*c^8 - 200*a^6*b*c^7 - 96*a^7*b*c^6 - 47*a^2*b^2*c^10 + 145*a^2
*b^3*c^9 - 354*a^2*b^4*c^8 + 612*a^2*b^5*c^7 - 655*a^2*b^6*c^6 + 635*a^2*b^7*c^5 - 202*a^2*b^8*c^4 - 222*a^2*b
^9*c^3 + 4*a^2*b^10*c^2 + 239*a^3*b^2*c^9 - 524*a^3*b^3*c^8 + 536*a^3*b^4*c^7 - 564*a^3*b^5*c^6 - 115*a^3*b^6*
c^5 + 856*a^3*b^7*c^4 + 2*a^3*b^8*c^3 - 20*a^3*b^9*c^2 - 37*a^4*b^2*c^8 + 9*a^4*b^3*c^7 + 583*a^4*b^4*c^6 - 13
62*a^4*b^5*c^5 - 152*a^4*b^6*c^4 + 156*a^4*b^7*c^3 + 8*a^4*b^8*c^2 - 399*a^5*b^2*c^7 + 904*a^5*b^3*c^6 + 394*a
^5*b^4*c^5 - 388*a^5*b^5*c^4 - 60*a^5*b^6*c^3 - 340*a^6*b^2*c^6 + 364*a^6*b^3*c^5 + 136*a^6*b^4*c^4 - 100*a^7*
b^2*c^5))/c^8 + (((2048*(16*a^3*c^13 - 32*a^2*c^14 + 176*a^4*c^12 + 176*a^5*c^11 + 48*a^6*c^10 - 2*b^4*c^12 +
6*b^5*c^11 - 18*b^6*c^10 + 26*b^7*c^9 - 12*b^8*c^8 + 16*a*b^2*c^13 - 40*a*b^3*c^12 + 122*a*b^4*c^11 - 192*a*b^
5*c^10 + 74*a*b^6*c^9 + 20*a*b^7*c^8 + 64*a^2*b*c^13 - 144*a^3*b*c^12 - 352*a^4*b*c^11 - 144*a^5*b*c^10 - 204*
a^2*b^2*c^12 + 388*a^2*b^3*c^11 - 50*a^2*b^4*c^10 - 182*a^2*b^5*c^9 + 4*a^2*b^6*c^8 - 260*a^3*b^2*c^11 + 496*a
^3*b^3*c^10 + 10*a^3*b^4*c^9 - 20*a^3*b^5*c^8 - 148*a^4*b^2*c^10 + 116*a^4*b^3*c^9 + 8*a^4*b^4*c^8 - 44*a^5*b^
2*c^9))/c^8 - (2048*tan(x/2)*((b^10 - a^2*b^8 - 8*a^5*c^5 - 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) + 10*a^3*
b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) + 52*a^2*b^6*c^2 - 96*a^3*b^4*c^3 + 66*a^4*b^2*c^4 - 33*a^4*b^4*c^2 +
38*a^5*b^2*c^3 - 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3
*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(
1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*
c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2)*(32*a*c^16 - 64*a^2*c^15 - 128*a^3*c^14 + 64*a^4*c^13 + 96*a^5*c^12
- 8*b^2*c^15 + 24*b^3*c^14 - 32*b^4*c^13 + 32*b^5*c^12 - 24*b^6*c^11 + 8*b^7*c^10 + 144*a*b^2*c^14 - 200*a*b^
3*c^13 + 184*a*b^4*c^12 - 56*a*b^5*c^11 - 8*a*b^6*c^10 + 288*a^2*b*c^14 + 352*a^3*b*c^13 - 32*a^4*b*c^12 - 320
*a^2*b^2*c^13 + 8*a^2*b^3*c^12 + 96*a^2*b^4*c^11 - 8*a^2*b^5*c^10 - 272*a^3*b^2*c^12 + 40*a^3*b^3*c^11 + 8*a^3
*b^4*c^10 - 56*a^4*b^2*c^11 - 96*a*b*c^15))/c^8)*((b^10 - a^2*b^8 - 8*a^5*c^5 - 8*a^6*c^4 - b^7*(-(4*a*c - b^2
)^3)^(1/2) + 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) + 52*a^2*b^6*c^2 - 96*a^3*b^4*c^3 + 66*a^4*b^2*c^
4 - 33*a^4*b^4*c^2 + 38*a^5*b^2*c^3 - 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c
- b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*
(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b
^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2) - (2048*tan(x/2)*(8*a*c^14 - 64*a^2*c^13 + 80*a
^3*c^12 + 168*a^4*c^11 - 192*a^5*c^10 - 136*a^6*c^9 + 72*a^7*c^8 - 2*b^2*c^13 + 6*b^3*c^12 - 17*b^4*c^11 + 33*
b^5*c^10 - 49*b^6*c^9 + 61*b^7*c^8 - 52*b^8*c^7 + 36*b^9*c^6 - 24*b^10*c^5 + 8*b^11*c^4 + 84*a*b^2*c^12 - 178*
a*b^3*c^11 + 295*a*b^4*c^10 - 416*a*b^5*c^9 + 375*a*b^6*c^8 - 308*a*b^7*c^7 + 244*a*b^8*c^6 - 72*a*b^9*c^5 - 8
*a*b^10*c^4 + 184*a^2*b*c^12 - 328*a^3*b*c^11 - 16*a^4*b*c^10 + 496*a^5*b*c^9 - 88*a^6*b*c^8 - 416*a^2*b^2*c^1
1 + 770*a^2*b^3*c^10 - 723*a^2*b^4*c^9 + 779*a^2*b^5*c^8 - 732*a^2*b^6*c^7 + 80*a^2*b^7*c^6 + 112*a^2*b^8*c^5
- 8*a^2*b^9*c^4 + 180*a^3*b^2*c^10 - 494*a^3*b^3*c^9 + 521*a^3*b^4*c^8 + 572*a^3*b^5*c^7 - 424*a^3*b^6*c^6 + 5
6*a^3*b^7*c^5 + 8*a^3*b^8*c^4 + 234*a^4*b^2*c^9 - 1152*a^4*b^3*c^8 + 416*a^4*b^4*c^7 - 140*a^4*b^5*c^6 - 72*a^
4*b^6*c^5 + 64*a^5*b^2*c^8 + 192*a^5*b^3*c^7 + 220*a^5*b^4*c^6 - 256*a^6*b^2*c^7 - 24*a*b*c^13))/c^8)*((b^10 -
a^2*b^8 - 8*a^5*c^5 - 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) + 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1
/2) + 52*a^2*b^6*c^2 - 96*a^3*b^4*c^3 + 66*a^4*b^2*c^4 - 33*a^4*b^4*c^2 + 38*a^5*b^2*c^3 - 12*a*b^8*c + 4*a^3*
b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) -
10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 +
16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^
(1/2))*((b^10 - a^2*b^8 - 8*a^5*c^5 - 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) + 10*a^3*b^6*c + a^2*b^5*(-(4*a
*c - b^2)^3)^(1/2) + 52*a^2*b^6*c^2 - 96*a^3*b^4*c^3 + 66*a^4*b^2*c^4 - 33*a^4*b^4*c^2 + 38*a^5*b^2*c^3 - 12*a
*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c -
b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10
+ 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*
a^3*b^2*c^7)))^(1/2) - (2048*tan(x/2)*(20*a*b^12 + 4*b^12*c - 4*b^13 - 40*a^2*b^11 + 40*a^3*b^10 - 20*a^4*b^9
+ 4*a^5*b^8 + 2*a^4*c^9 - 18*a^5*c^8 + 38*a^6*c^7 + 2*a^7*c^6 - 44*a^8*c^5 + 12*a^9*c^4 + b^8*c^5 - b^9*c^4 +
4*b^10*c^3 - 4*b^11*c^2 - 8*a*b^6*c^6 + 4*a*b^7*c^5 - 31*a*b^8*c^4 + 20*a*b^9*c^3 - 20*a*b^10*c^2 - 160*a^2*b^
10*c + 320*a^3*b^9*c + 26*a^4*b*c^8 - 300*a^4*b^8*c - 84*a^5*b*c^7 + 136*a^5*b^7*c + 2*a^6*b*c^6 - 24*a^6*b^6*
c + 168*a^7*b*c^5 - 92*a^8*b*c^4 + 20*a^2*b^4*c^7 + 8*a^2*b^5*c^6 + 82*a^2*b^6*c^5 + 6*a^2*b^7*c^4 + 8*a^2*b^8
*c^3 - 44*a^2*b^9*c^2 - 16*a^3*b^2*c^8 - 40*a^3*b^3*c^7 - 104*a^3*b^4*c^6 - 132*a^3*b^5*c^5 + 34*a^3*b^6*c^4 +
72*a^3*b^7*c^3 + 460*a^3*b^8*c^2 + 82*a^4*b^2*c^7 + 174*a^4*b^3*c^6 + 41*a^4*b^4*c^5 - 149*a^4*b^5*c^4 - 660*
a^4*b^6*c^3 - 900*a^4*b^7*c^2 - 90*a^5*b^2*c^6 + 96*a^5*b^3*c^5 + 541*a^5*b^4*c^4 + 1156*a^5*b^5*c^3 + 764*a^5
*b^6*c^2 - 204*a^6*b^2*c^5 - 704*a^6*b^3*c^4 - 840*a^6*b^4*c^3 - 300*a^6*b^5*c^2 + 384*a^7*b^2*c^4 + 272*a^7*b
^3*c^3 + 44*a^7*b^4*c^2 - 32*a^8*b^2*c^3 + 24*a*b^11*c))/c^8)*((b^10 - a^2*b^8 - 8*a^5*c^5 - 8*a^6*c^4 - b^7*(
-(4*a*c - b^2)^3)^(1/2) + 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) + 52*a^2*b^6*c^2 - 96*a^3*b^4*c^3 +
66*a^4*b^2*c^4 - 33*a^4*b^4*c^2 + 38*a^5*b^2*c^3 - 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b
^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2)
+ 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2
*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2) + (((2048*(12*a^3*c^11 - 28*a^4*c^
10 - 44*a^5*c^9 + 72*a^6*c^8 + 88*a^7*c^7 + 12*a^8*c^6 + b^5*c^9 - 4*b^6*c^8 + 10*b^7*c^7 - 20*b^8*c^6 + 29*b^
9*c^5 - 30*b^10*c^4 + 26*b^11*c^3 - 12*b^12*c^2 - 6*a*b^3*c^10 + 27*a*b^4*c^9 - 72*a*b^5*c^8 + 154*a*b^6*c^7 -
238*a*b^7*c^6 + 251*a*b^8*c^5 - 228*a*b^9*c^4 + 98*a*b^10*c^3 + 20*a*b^11*c^2 + 8*a^2*b*c^11 - 68*a^3*b*c^10
+ 112*a^4*b*c^9 + 100*a^5*b*c^8 - 200*a^6*b*c^7 - 96*a^7*b*c^6 - 47*a^2*b^2*c^10 + 145*a^2*b^3*c^9 - 354*a^2*b
^4*c^8 + 612*a^2*b^5*c^7 - 655*a^2*b^6*c^6 + 635*a^2*b^7*c^5 - 202*a^2*b^8*c^4 - 222*a^2*b^9*c^3 + 4*a^2*b^10*
c^2 + 239*a^3*b^2*c^9 - 524*a^3*b^3*c^8 + 536*a^3*b^4*c^7 - 564*a^3*b^5*c^6 - 115*a^3*b^6*c^5 + 856*a^3*b^7*c^
4 + 2*a^3*b^8*c^3 - 20*a^3*b^9*c^2 - 37*a^4*b^2*c^8 + 9*a^4*b^3*c^7 + 583*a^4*b^4*c^6 - 1362*a^4*b^5*c^5 - 152
*a^4*b^6*c^4 + 156*a^4*b^7*c^3 + 8*a^4*b^8*c^2 - 399*a^5*b^2*c^7 + 904*a^5*b^3*c^6 + 394*a^5*b^4*c^5 - 388*a^5
*b^5*c^4 - 60*a^5*b^6*c^3 - 340*a^6*b^2*c^6 + 364*a^6*b^3*c^5 + 136*a^6*b^4*c^4 - 100*a^7*b^2*c^5))/c^8 + (((2
048*(16*a^3*c^13 - 32*a^2*c^14 + 176*a^4*c^12 + 176*a^5*c^11 + 48*a^6*c^10 - 2*b^4*c^12 + 6*b^5*c^11 - 18*b^6*
c^10 + 26*b^7*c^9 - 12*b^8*c^8 + 16*a*b^2*c^13 - 40*a*b^3*c^12 + 122*a*b^4*c^11 - 192*a*b^5*c^10 + 74*a*b^6*c^
9 + 20*a*b^7*c^8 + 64*a^2*b*c^13 - 144*a^3*b*c^12 - 352*a^4*b*c^11 - 144*a^5*b*c^10 - 204*a^2*b^2*c^12 + 388*a
^2*b^3*c^11 - 50*a^2*b^4*c^10 - 182*a^2*b^5*c^9 + 4*a^2*b^6*c^8 - 260*a^3*b^2*c^11 + 496*a^3*b^3*c^10 + 10*a^3
*b^4*c^9 - 20*a^3*b^5*c^8 - 148*a^4*b^2*c^10 + 116*a^4*b^3*c^9 + 8*a^4*b^4*c^8 - 44*a^5*b^2*c^9))/c^8 + (2048*
tan(x/2)*((b^10 - a^2*b^8 - 8*a^5*c^5 - 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) + 10*a^3*b^6*c + a^2*b^5*(-(4
*a*c - b^2)^3)^(1/2) + 52*a^2*b^6*c^2 - 96*a^3*b^4*c^3 + 66*a^4*b^2*c^4 - 33*a^4*b^4*c^2 + 38*a^5*b^2*c^3 - 12
*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c
- b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^1
0 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 -
8*a^3*b^2*c^7)))^(1/2)*(32*a*c^16 - 64*a^2*c^15 - 128*a^3*c^14 + 64*a^4*c^13 + 96*a^5*c^12 - 8*b^2*c^15 + 24*b
^3*c^14 - 32*b^4*c^13 + 32*b^5*c^12 - 24*b^6*c^11 + 8*b^7*c^10 + 144*a*b^2*c^14 - 200*a*b^3*c^13 + 184*a*b^4*c
^12 - 56*a*b^5*c^11 - 8*a*b^6*c^10 + 288*a^2*b*c^14 + 352*a^3*b*c^13 - 32*a^4*b*c^12 - 320*a^2*b^2*c^13 + 8*a^
2*b^3*c^12 + 96*a^2*b^4*c^11 - 8*a^2*b^5*c^10 - 272*a^3*b^2*c^12 + 40*a^3*b^3*c^11 + 8*a^3*b^4*c^10 - 56*a^4*b
^2*c^11 - 96*a*b*c^15))/c^8)*((b^10 - a^2*b^8 - 8*a^5*c^5 - 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) + 10*a^3*
b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) + 52*a^2*b^6*c^2 - 96*a^3*b^4*c^3 + 66*a^4*b^2*c^4 - 33*a^4*b^4*c^2 +
38*a^5*b^2*c^3 - 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3
*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(
1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*
c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2) + (2048*tan(x/2)*(8*a*c^14 - 64*a^2*c^13 + 80*a^3*c^12 + 168*a^4*c^
11 - 192*a^5*c^10 - 136*a^6*c^9 + 72*a^7*c^8 - 2*b^2*c^13 + 6*b^3*c^12 - 17*b^4*c^11 + 33*b^5*c^10 - 49*b^6*c^
9 + 61*b^7*c^8 - 52*b^8*c^7 + 36*b^9*c^6 - 24*b^10*c^5 + 8*b^11*c^4 + 84*a*b^2*c^12 - 178*a*b^3*c^11 + 295*a*b
^4*c^10 - 416*a*b^5*c^9 + 375*a*b^6*c^8 - 308*a*b^7*c^7 + 244*a*b^8*c^6 - 72*a*b^9*c^5 - 8*a*b^10*c^4 + 184*a^
2*b*c^12 - 328*a^3*b*c^11 - 16*a^4*b*c^10 + 496*a^5*b*c^9 - 88*a^6*b*c^8 - 416*a^2*b^2*c^11 + 770*a^2*b^3*c^10
- 723*a^2*b^4*c^9 + 779*a^2*b^5*c^8 - 732*a^2*b^6*c^7 + 80*a^2*b^7*c^6 + 112*a^2*b^8*c^5 - 8*a^2*b^9*c^4 + 18
0*a^3*b^2*c^10 - 494*a^3*b^3*c^9 + 521*a^3*b^4*c^8 + 572*a^3*b^5*c^7 - 424*a^3*b^6*c^6 + 56*a^3*b^7*c^5 + 8*a^
3*b^8*c^4 + 234*a^4*b^2*c^9 - 1152*a^4*b^3*c^8 + 416*a^4*b^4*c^7 - 140*a^4*b^5*c^6 - 72*a^4*b^6*c^5 + 64*a^5*b
^2*c^8 + 192*a^5*b^3*c^7 + 220*a^5*b^4*c^6 - 256*a^6*b^2*c^7 - 24*a*b*c^13))/c^8)*((b^10 - a^2*b^8 - 8*a^5*c^5
- 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) + 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) + 52*a^2*b^6*c^2
- 96*a^3*b^4*c^3 + 66*a^4*b^2*c^4 - 33*a^4*b^4*c^2 + 38*a^5*b^2*c^3 - 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2
)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4
*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8
- b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2))*((b^10 - a^2*
b^8 - 8*a^5*c^5 - 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1/2) + 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) +
52*a^2*b^6*c^2 - 96*a^3*b^4*c^3 + 66*a^4*b^2*c^4 - 33*a^4*b^4*c^2 + 38*a^5*b^2*c^3 - 12*a*b^8*c + 4*a^3*b*c^3
*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a
^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^
4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2)
+ (2048*tan(x/2)*(20*a*b^12 + 4*b^12*c - 4*b^13 - 40*a^2*b^11 + 40*a^3*b^10 - 20*a^4*b^9 + 4*a^5*b^8 + 2*a^4*
c^9 - 18*a^5*c^8 + 38*a^6*c^7 + 2*a^7*c^6 - 44*a^8*c^5 + 12*a^9*c^4 + b^8*c^5 - b^9*c^4 + 4*b^10*c^3 - 4*b^11*
c^2 - 8*a*b^6*c^6 + 4*a*b^7*c^5 - 31*a*b^8*c^4 + 20*a*b^9*c^3 - 20*a*b^10*c^2 - 160*a^2*b^10*c + 320*a^3*b^9*c
+ 26*a^4*b*c^8 - 300*a^4*b^8*c - 84*a^5*b*c^7 + 136*a^5*b^7*c + 2*a^6*b*c^6 - 24*a^6*b^6*c + 168*a^7*b*c^5 -
92*a^8*b*c^4 + 20*a^2*b^4*c^7 + 8*a^2*b^5*c^6 + 82*a^2*b^6*c^5 + 6*a^2*b^7*c^4 + 8*a^2*b^8*c^3 - 44*a^2*b^9*c^
2 - 16*a^3*b^2*c^8 - 40*a^3*b^3*c^7 - 104*a^3*b^4*c^6 - 132*a^3*b^5*c^5 + 34*a^3*b^6*c^4 + 72*a^3*b^7*c^3 + 46
0*a^3*b^8*c^2 + 82*a^4*b^2*c^7 + 174*a^4*b^3*c^6 + 41*a^4*b^4*c^5 - 149*a^4*b^5*c^4 - 660*a^4*b^6*c^3 - 900*a^
4*b^7*c^2 - 90*a^5*b^2*c^6 + 96*a^5*b^3*c^5 + 541*a^5*b^4*c^4 + 1156*a^5*b^5*c^3 + 764*a^5*b^6*c^2 - 204*a^6*b
^2*c^5 - 704*a^6*b^3*c^4 - 840*a^6*b^4*c^3 - 300*a^6*b^5*c^2 + 384*a^7*b^2*c^4 + 272*a^7*b^3*c^3 + 44*a^7*b^4*
c^2 - 32*a^8*b^2*c^3 + 24*a*b^11*c))/c^8)*((b^10 - a^2*b^8 - 8*a^5*c^5 - 8*a^6*c^4 - b^7*(-(4*a*c - b^2)^3)^(1
/2) + 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) + 52*a^2*b^6*c^2 - 96*a^3*b^4*c^3 + 66*a^4*b^2*c^4 - 33*
a^4*b^4*c^2 + 38*a^5*b^2*c^3 - 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c*(-(4*a*c - b^2)
^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6*a*b^5*c*(-(4*a*
c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9 + 10*a*b^4*c^7
- 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2)))*((b^10 - a^2*b^8 - 8*a^5*c^5 - 8*a^6*c^4 - b^7*(-(4*
a*c - b^2)^3)^(1/2) + 10*a^3*b^6*c + a^2*b^5*(-(4*a*c - b^2)^3)^(1/2) + 52*a^2*b^6*c^2 - 96*a^3*b^4*c^3 + 66*a
^4*b^2*c^4 - 33*a^4*b^4*c^2 + 38*a^5*b^2*c^3 - 12*a*b^8*c + 4*a^3*b*c^3*(-(4*a*c - b^2)^3)^(1/2) - 4*a^3*b^3*c
*(-(4*a*c - b^2)^3)^(1/2) + 3*a^4*b*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a^2*b^3*c^2*(-(4*a*c - b^2)^3)^(1/2) + 6
*a*b^5*c*(-(4*a*c - b^2)^3)^(1/2))/(2*(16*a^2*c^10 + 32*a^3*c^9 + 16*a^4*c^8 + b^4*c^8 - b^6*c^6 - 8*a*b^2*c^9
+ 10*a*b^4*c^7 - 32*a^2*b^2*c^8 + a^2*b^4*c^6 - 8*a^3*b^2*c^7)))^(1/2)*2i + (atan(-((((2048*tan(x/2)*(20*a*b^
12 + 4*b^12*c - 4*b^13 - 40*a^2*b^11 + 40*a^3*b^10 - 20*a^4*b^9 + 4*a^5*b^8 + 2*a^4*c^9 - 18*a^5*c^8 + 38*a^6*
c^7 + 2*a^7*c^6 - 44*a^8*c^5 + 12*a^9*c^4 + b^8*c^5 - b^9*c^4 + 4*b^10*c^3 - 4*b^11*c^2 - 8*a*b^6*c^6 + 4*a*b^
7*c^5 - 31*a*b^8*c^4 + 20*a*b^9*c^3 - 20*a*b^10*c^2 - 160*a^2*b^10*c + 320*a^3*b^9*c + 26*a^4*b*c^8 - 300*a^4*
b^8*c - 84*a^5*b*c^7 + 136*a^5*b^7*c + 2*a^6*b*c^6 - 24*a^6*b^6*c + 168*a^7*b*c^5 - 92*a^8*b*c^4 + 20*a^2*b^4*
c^7 + 8*a^2*b^5*c^6 + 82*a^2*b^6*c^5 + 6*a^2*b^7*c^4 + 8*a^2*b^8*c^3 - 44*a^2*b^9*c^2 - 16*a^3*b^2*c^8 - 40*a^
3*b^3*c^7 - 104*a^3*b^4*c^6 - 132*a^3*b^5*c^5 + 34*a^3*b^6*c^4 + 72*a^3*b^7*c^3 + 460*a^3*b^8*c^2 + 82*a^4*b^2
*c^7 + 174*a^4*b^3*c^6 + 41*a^4*b^4*c^5 - 149*a^4*b^5*c^4 - 660*a^4*b^6*c^3 - 900*a^4*b^7*c^2 - 90*a^5*b^2*c^6
+ 96*a^5*b^3*c^5 + 541*a^5*b^4*c^4 + 1156*a^5*b^5*c^3 + 764*a^5*b^6*c^2 - 204*a^6*b^2*c^5 - 704*a^6*b^3*c^4 -
840*a^6*b^4*c^3 - 300*a^6*b^5*c^2 + 384*a^7*b^2*c^4 + 272*a^7*b^3*c^3 + 44*a^7*b^4*c^2 - 32*a^8*b^2*c^3 + 24*
a*b^11*c))/c^8 - (((2048*(12*a^3*c^11 - 28*a^4*c^10 - 44*a^5*c^9 + 72*a^6*c^8 + 88*a^7*c^7 + 12*a^8*c^6 + b^5*
c^9 - 4*b^6*c^8 + 10*b^7*c^7 - 20*b^8*c^6 + 29*b^9*c^5 - 30*b^10*c^4 + 26*b^11*c^3 - 12*b^12*c^2 - 6*a*b^3*c^1
0 + 27*a*b^4*c^9 - 72*a*b^5*c^8 + 154*a*b^6*c^7 - 238*a*b^7*c^6 + 251*a*b^8*c^5 - 228*a*b^9*c^4 + 98*a*b^10*c^
3 + 20*a*b^11*c^2 + 8*a^2*b*c^11 - 68*a^3*b*c^10 + 112*a^4*b*c^9 + 100*a^5*b*c^8 - 200*a^6*b*c^7 - 96*a^7*b*c^
6 - 47*a^2*b^2*c^10 + 145*a^2*b^3*c^9 - 354*a^2*b^4*c^8 + 612*a^2*b^5*c^7 - 655*a^2*b^6*c^6 + 635*a^2*b^7*c^5
- 202*a^2*b^8*c^4 - 222*a^2*b^9*c^3 + 4*a^2*b^10*c^2 + 239*a^3*b^2*c^9 - 524*a^3*b^3*c^8 + 536*a^3*b^4*c^7 - 5
64*a^3*b^5*c^6 - 115*a^3*b^6*c^5 + 856*a^3*b^7*c^4 + 2*a^3*b^8*c^3 - 20*a^3*b^9*c^2 - 37*a^4*b^2*c^8 + 9*a^4*b
^3*c^7 + 583*a^4*b^4*c^6 - 1362*a^4*b^5*c^5 - 152*a^4*b^6*c^4 + 156*a^4*b^7*c^3 + 8*a^4*b^8*c^2 - 399*a^5*b^2*
c^7 + 904*a^5*b^3*c^6 + 394*a^5*b^4*c^5 - 388*a^5*b^5*c^4 - 60*a^5*b^6*c^3 - 340*a^6*b^2*c^6 + 364*a^6*b^3*c^5
+ 136*a^6*b^4*c^4 - 100*a^7*b^2*c^5))/c^8 - (((2048*tan(x/2)*(8*a*c^14 - 64*a^2*c^13 + 80*a^3*c^12 + 168*a^4*
c^11 - 192*a^5*c^10 - 136*a^6*c^9 + 72*a^7*c^8 - 2*b^2*c^13 + 6*b^3*c^12 - 17*b^4*c^11 + 33*b^5*c^10 - 49*b^6*
c^9 + 61*b^7*c^8 - 52*b^8*c^7 + 36*b^9*c^6 - 24*b^10*c^5 + 8*b^11*c^4 + 84*a*b^2*c^12 - 178*a*b^3*c^11 + 295*a
*b^4*c^10 - 416*a*b^5*c^9 + 375*a*b^6*c^8 - 308*a*b^7*c^7 + 244*a*b^8*c^6 - 72*a*b^9*c^5 - 8*a*b^10*c^4 + 184*
a^2*b*c^12 - 328*a^3*b*c^11 - 16*a^4*b*c^10 + 496*a^5*b*c^9 - 88*a^6*b*c^8 - 416*a^2*b^2*c^11 + 770*a^2*b^3*c^
10 - 723*a^2*b^4*c^9 + 779*a^2*b^5*c^8 - 732*a^2*b^6*c^7 + 80*a^2*b^7*c^6 + 112*a^2*b^8*c^5 - 8*a^2*b^9*c^4 +
180*a^3*b^2*c^10 - 494*a^3*b^3*c^9 + 521*a^3*b^4*c^8 + 572*a^3*b^5*c^7 - 424*a^3*b^6*c^6 + 56*a^3*b^7*c^5 + 8*
a^3*b^8*c^4 + 234*a^4*b^2*c^9 - 1152*a^4*b^3*c^8 + 416*a^4*b^4*c^7 - 140*a^4*b^5*c^6 - 72*a^4*b^6*c^5 + 64*a^5
*b^2*c^8 + 192*a^5*b^3*c^7 + 220*a^5*b^4*c^6 - 256*a^6*b^2*c^7 - 24*a*b*c^13))/c^8 - (((2048*(16*a^3*c^13 - 32
*a^2*c^14 + 176*a^4*c^12 + 176*a^5*c^11 + 48*a^6*c^10 - 2*b^4*c^12 + 6*b^5*c^11 - 18*b^6*c^10 + 26*b^7*c^9 - 1
2*b^8*c^8 + 16*a*b^2*c^13 - 40*a*b^3*c^12 + 122*a*b^4*c^11 - 192*a*b^5*c^10 + 74*a*b^6*c^9 + 20*a*b^7*c^8 + 64
*a^2*b*c^13 - 144*a^3*b*c^12 - 352*a^4*b*c^11 - 144*a^5*b*c^10 - 204*a^2*b^2*c^12 + 388*a^2*b^3*c^11 - 50*a^2*
b^4*c^10 - 182*a^2*b^5*c^9 + 4*a^2*b^6*c^8 - 260*a^3*b^2*c^11 + 496*a^3*b^3*c^10 + 10*a^3*b^4*c^9 - 20*a^3*b^5
*c^8 - 148*a^4*b^2*c^10 + 116*a^4*b^3*c^9 + 8*a^4*b^4*c^8 - 44*a^5*b^2*c^9))/c^8 - (1024*tan(x/2)*(b^2*2i - a*
c*2i + c^2*1i)*(32*a*c^16 - 64*a^2*c^15 - 128*a^3*c^14 + 64*a^4*c^13 + 96*a^5*c^12 - 8*b^2*c^15 + 24*b^3*c^14
- 32*b^4*c^13 + 32*b^5*c^12 - 24*b^6*c^11 + 8*b^7*c^10 + 144*a*b^2*c^14 - 200*a*b^3*c^13 + 184*a*b^4*c^12 - 56
*a*b^5*c^11 - 8*a*b^6*c^10 + 288*a^2*b*c^14 + 352*a^3*b*c^13 - 32*a^4*b*c^12 - 320*a^2*b^2*c^13 + 8*a^2*b^3*c^
12 + 96*a^2*b^4*c^11 - 8*a^2*b^5*c^10 - 272*a^3*b^2*c^12 + 40*a^3*b^3*c^11 + 8*a^3*b^4*c^10 - 56*a^4*b^2*c^11
- 96*a*b*c^15))/c^11)*(b^2*2i - a*c*2i + c^2*1i))/(2*c^3))*(b^2*2i - a*c*2i + c^2*1i))/(2*c^3))*(b^2*2i - a*c*
2i + c^2*1i))/(2*c^3))*(b^2*2i - a*c*2i + c^2*1i)*1i)/(2*c^3) + (((2048*tan(x/2)*(20*a*b^12 + 4*b^12*c - 4*b^1
3 - 40*a^2*b^11 + 40*a^3*b^10 - 20*a^4*b^9 + 4*a^5*b^8 + 2*a^4*c^9 - 18*a^5*c^8 + 38*a^6*c^7 + 2*a^7*c^6 - 44*
a^8*c^5 + 12*a^9*c^4 + b^8*c^5 - b^9*c^4 + 4*b^10*c^3 - 4*b^11*c^2 - 8*a*b^6*c^6 + 4*a*b^7*c^5 - 31*a*b^8*c^4
+ 20*a*b^9*c^3 - 20*a*b^10*c^2 - 160*a^2*b^10*c + 320*a^3*b^9*c + 26*a^4*b*c^8 - 300*a^4*b^8*c - 84*a^5*b*c^7
+ 136*a^5*b^7*c + 2*a^6*b*c^6 - 24*a^6*b^6*c + 168*a^7*b*c^5 - 92*a^8*b*c^4 + 20*a^2*b^4*c^7 + 8*a^2*b^5*c^6 +
82*a^2*b^6*c^5 + 6*a^2*b^7*c^4 + 8*a^2*b^8*c^3 - 44*a^2*b^9*c^2 - 16*a^3*b^2*c^8 - 40*a^3*b^3*c^7 - 104*a^3*b
^4*c^6 - 132*a^3*b^5*c^5 + 34*a^3*b^6*c^4 + 72*a^3*b^7*c^3 + 460*a^3*b^8*c^2 + 82*a^4*b^2*c^7 + 174*a^4*b^3*c^
6 + 41*a^4*b^4*c^5 - 149*a^4*b^5*c^4 - 660*a^4*b^6*c^3 - 900*a^4*b^7*c^2 - 90*a^5*b^2*c^6 + 96*a^5*b^3*c^5 + 5
41*a^5*b^4*c^4 + 1156*a^5*b^5*c^3 + 764*a^5*b^6*c^2 - 204*a^6*b^2*c^5 - 704*a^6*b^3*c^4 - 840*a^6*b^4*c^3 - 30
0*a^6*b^5*c^2 + 384*a^7*b^2*c^4 + 272*a^7*b^3*c^3 + 44*a^7*b^4*c^2 - 32*a^8*b^2*c^3 + 24*a*b^11*c))/c^8 + (((2
048*(12*a^3*c^11 - 28*a^4*c^10 - 44*a^5*c^9 + 72*a^6*c^8 + 88*a^7*c^7 + 12*a^8*c^6 + b^5*c^9 - 4*b^6*c^8 + 10*
b^7*c^7 - 20*b^8*c^6 + 29*b^9*c^5 - 30*b^10*c^4 + 26*b^11*c^3 - 12*b^12*c^2 - 6*a*b^3*c^10 + 27*a*b^4*c^9 - 72
*a*b^5*c^8 + 154*a*b^6*c^7 - 238*a*b^7*c^6 + 251*a*b^8*c^5 - 228*a*b^9*c^4 + 98*a*b^10*c^3 + 20*a*b^11*c^2 + 8
*a^2*b*c^11 - 68*a^3*b*c^10 + 112*a^4*b*c^9 + 100*a^5*b*c^8 - 200*a^6*b*c^7 - 96*a^7*b*c^6 - 47*a^2*b^2*c^10 +
145*a^2*b^3*c^9 - 354*a^2*b^4*c^8 + 612*a^2*b^5*c^7 - 655*a^2*b^6*c^6 + 635*a^2*b^7*c^5 - 202*a^2*b^8*c^4 - 2
22*a^2*b^9*c^3 + 4*a^2*b^10*c^2 + 239*a^3*b^2*c^9 - 524*a^3*b^3*c^8 + 536*a^3*b^4*c^7 - 564*a^3*b^5*c^6 - 115*
a^3*b^6*c^5 + 856*a^3*b^7*c^4 + 2*a^3*b^8*c^3 - 20*a^3*b^9*c^2 - 37*a^4*b^2*c^8 + 9*a^4*b^3*c^7 + 583*a^4*b^4*
c^6 - 1362*a^4*b^5*c^5 - 152*a^4*b^6*c^4 + 156*a^4*b^7*c^3 + 8*a^4*b^8*c^2 - 399*a^5*b^2*c^7 + 904*a^5*b^3*c^6
+ 394*a^5*b^4*c^5 - 388*a^5*b^5*c^4 - 60*a^5*b^6*c^3 - 340*a^6*b^2*c^6 + 364*a^6*b^3*c^5 + 136*a^6*b^4*c^4 -
100*a^7*b^2*c^5))/c^8 + (((2048*tan(x/2)*(8*a*c^14 - 64*a^2*c^13 + 80*a^3*c^12 + 168*a^4*c^11 - 192*a^5*c^10 -
136*a^6*c^9 + 72*a^7*c^8 - 2*b^2*c^13 + 6*b^3*c^12 - 17*b^4*c^11 + 33*b^5*c^10 - 49*b^6*c^9 + 61*b^7*c^8 - 52
*b^8*c^7 + 36*b^9*c^6 - 24*b^10*c^5 + 8*b^11*c^4 + 84*a*b^2*c^12 - 178*a*b^3*c^11 + 295*a*b^4*c^10 - 416*a*b^5
*c^9 + 375*a*b^6*c^8 - 308*a*b^7*c^7 + 244*a*b^8*c^6 - 72*a*b^9*c^5 - 8*a*b^10*c^4 + 184*a^2*b*c^12 - 328*a^3*
b*c^11 - 16*a^4*b*c^10 + 496*a^5*b*c^9 - 88*a^6*b*c^8 - 416*a^2*b^2*c^11 + 770*a^2*b^3*c^10 - 723*a^2*b^4*c^9
+ 779*a^2*b^5*c^8 - 732*a^2*b^6*c^7 + 80*a^2*b^7*c^6 + 112*a^2*b^8*c^5 - 8*a^2*b^9*c^4 + 180*a^3*b^2*c^10 - 49
4*a^3*b^3*c^9 + 521*a^3*b^4*c^8 + 572*a^3*b^5*c^7 - 424*a^3*b^6*c^6 + 56*a^3*b^7*c^5 + 8*a^3*b^8*c^4 + 234*a^4
*b^2*c^9 - 1152*a^4*b^3*c^8 + 416*a^4*b^4*c^7 - 140*a^4*b^5*c^6 - 72*a^4*b^6*c^5 + 64*a^5*b^2*c^8 + 192*a^5*b^
3*c^7 + 220*a^5*b^4*c^6 - 256*a^6*b^2*c^7 - 24*a*b*c^13))/c^8 + (((2048*(16*a^3*c^13 - 32*a^2*c^14 + 176*a^4*c
^12 + 176*a^5*c^11 + 48*a^6*c^10 - 2*b^4*c^12 + 6*b^5*c^11 - 18*b^6*c^10 + 26*b^7*c^9 - 12*b^8*c^8 + 16*a*b^2*
c^13 - 40*a*b^3*c^12 + 122*a*b^4*c^11 - 192*a*b^5*c^10 + 74*a*b^6*c^9 + 20*a*b^7*c^8 + 64*a^2*b*c^13 - 144*a^3
*b*c^12 - 352*a^4*b*c^11 - 144*a^5*b*c^10 - 204*a^2*b^2*c^12 + 388*a^2*b^3*c^11 - 50*a^2*b^4*c^10 - 182*a^2*b^
5*c^9 + 4*a^2*b^6*c^8 - 260*a^3*b^2*c^11 + 496*a^3*b^3*c^10 + 10*a^3*b^4*c^9 - 20*a^3*b^5*c^8 - 148*a^4*b^2*c^
10 + 116*a^4*b^3*c^9 + 8*a^4*b^4*c^8 - 44*a^5*b^2*c^9))/c^8 + (1024*tan(x/2)*(b^2*2i - a*c*2i + c^2*1i)*(32*a*
c^16 - 64*a^2*c^15 - 128*a^3*c^14 + 64*a^4*c^13 + 96*a^5*c^12 - 8*b^2*c^15 + 24*b^3*c^14 - 32*b^4*c^13 + 32*b^
5*c^12 - 24*b^6*c^11 + 8*b^7*c^10 + 144*a*b^2*c^14 - 200*a*b^3*c^13 + 184*a*b^4*c^12 - 56*a*b^5*c^11 - 8*a*b^6
*c^10 + 288*a^2*b*c^14 + 352*a^3*b*c^13 - 32*a^4*b*c^12 - 320*a^2*b^2*c^13 + 8*a^2*b^3*c^12 + 96*a^2*b^4*c^11
- 8*a^2*b^5*c^10 - 272*a^3*b^2*c^12 + 40*a^3*b^3*c^11 + 8*a^3*b^4*c^10 - 56*a^4*b^2*c^11 - 96*a*b*c^15))/c^11)
*(b^2*2i - a*c*2i + c^2*1i))/(2*c^3))*(b^2*2i - a*c*2i + c^2*1i))/(2*c^3))*(b^2*2i - a*c*2i + c^2*1i))/(2*c^3)
)*(b^2*2i - a*c*2i + c^2*1i)*1i)/(2*c^3))/((4096*(16*a^5*b^7 - 4*a^4*b^8 - 24*a^6*b^6 + 16*a^7*b^5 - 4*a^8*b^4
+ 3*a^6*c^6 - 10*a^7*c^5 + a^8*c^4 + 14*a^9*c^3 + 4*a^4*b^7*c - 2*a^5*b*c^6 + 4*a^5*b^6*c + 6*a^6*b*c^5 - 40*
a^6*b^5*c + 4*a^7*b*c^4 + 56*a^7*b^4*c - 22*a^8*b*c^3 - 28*a^8*b^3*c + 12*a^9*b*c^2 + 4*a^9*b^2*c + a^4*b^3*c^
5 - a^4*b^4*c^4 + 4*a^4*b^5*c^3 - 4*a^4*b^6*c^2 - a^5*b^2*c^5 - 8*a^5*b^3*c^4 + 10*a^6*b^2*c^4 - 4*a^6*b^3*c^3
- 8*a^6*b^4*c^2 + 4*a^7*b^2*c^3 + 48*a^7*b^3*c^2 - 48*a^8*b^2*c^2))/c^8 - (((2048*tan(x/2)*(20*a*b^12 + 4*b^1
2*c - 4*b^13 - 40*a^2*b^11 + 40*a^3*b^10 - 20*a^4*b^9 + 4*a^5*b^8 + 2*a^4*c^9 - 18*a^5*c^8 + 38*a^6*c^7 + 2*a^
7*c^6 - 44*a^8*c^5 + 12*a^9*c^4 + b^8*c^5 - b^9*c^4 + 4*b^10*c^3 - 4*b^11*c^2 - 8*a*b^6*c^6 + 4*a*b^7*c^5 - 31
*a*b^8*c^4 + 20*a*b^9*c^3 - 20*a*b^10*c^2 - 160*a^2*b^10*c + 320*a^3*b^9*c + 26*a^4*b*c^8 - 300*a^4*b^8*c - 84
*a^5*b*c^7 + 136*a^5*b^7*c + 2*a^6*b*c^6 - 24*a^6*b^6*c + 168*a^7*b*c^5 - 92*a^8*b*c^4 + 20*a^2*b^4*c^7 + 8*a^
2*b^5*c^6 + 82*a^2*b^6*c^5 + 6*a^2*b^7*c^4 + 8*a^2*b^8*c^3 - 44*a^2*b^9*c^2 - 16*a^3*b^2*c^8 - 40*a^3*b^3*c^7
- 104*a^3*b^4*c^6 - 132*a^3*b^5*c^5 + 34*a^3*b^6*c^4 + 72*a^3*b^7*c^3 + 460*a^3*b^8*c^2 + 82*a^4*b^2*c^7 + 174
*a^4*b^3*c^6 + 41*a^4*b^4*c^5 - 149*a^4*b^5*c^4 - 660*a^4*b^6*c^3 - 900*a^4*b^7*c^2 - 90*a^5*b^2*c^6 + 96*a^5*
b^3*c^5 + 541*a^5*b^4*c^4 + 1156*a^5*b^5*c^3 + 764*a^5*b^6*c^2 - 204*a^6*b^2*c^5 - 704*a^6*b^3*c^4 - 840*a^6*b
^4*c^3 - 300*a^6*b^5*c^2 + 384*a^7*b^2*c^4 + 272*a^7*b^3*c^3 + 44*a^7*b^4*c^2 - 32*a^8*b^2*c^3 + 24*a*b^11*c))
/c^8 - (((2048*(12*a^3*c^11 - 28*a^4*c^10 - 44*a^5*c^9 + 72*a^6*c^8 + 88*a^7*c^7 + 12*a^8*c^6 + b^5*c^9 - 4*b^
6*c^8 + 10*b^7*c^7 - 20*b^8*c^6 + 29*b^9*c^5 - 30*b^10*c^4 + 26*b^11*c^3 - 12*b^12*c^2 - 6*a*b^3*c^10 + 27*a*b
^4*c^9 - 72*a*b^5*c^8 + 154*a*b^6*c^7 - 238*a*b^7*c^6 + 251*a*b^8*c^5 - 228*a*b^9*c^4 + 98*a*b^10*c^3 + 20*a*b
^11*c^2 + 8*a^2*b*c^11 - 68*a^3*b*c^10 + 112*a^4*b*c^9 + 100*a^5*b*c^8 - 200*a^6*b*c^7 - 96*a^7*b*c^6 - 47*a^2
*b^2*c^10 + 145*a^2*b^3*c^9 - 354*a^2*b^4*c^8 + 612*a^2*b^5*c^7 - 655*a^2*b^6*c^6 + 635*a^2*b^7*c^5 - 202*a^2*
b^8*c^4 - 222*a^2*b^9*c^3 + 4*a^2*b^10*c^2 + 239*a^3*b^2*c^9 - 524*a^3*b^3*c^8 + 536*a^3*b^4*c^7 - 564*a^3*b^5
*c^6 - 115*a^3*b^6*c^5 + 856*a^3*b^7*c^4 + 2*a^3*b^8*c^3 - 20*a^3*b^9*c^2 - 37*a^4*b^2*c^8 + 9*a^4*b^3*c^7 + 5
83*a^4*b^4*c^6 - 1362*a^4*b^5*c^5 - 152*a^4*b^6*c^4 + 156*a^4*b^7*c^3 + 8*a^4*b^8*c^2 - 399*a^5*b^2*c^7 + 904*
a^5*b^3*c^6 + 394*a^5*b^4*c^5 - 388*a^5*b^5*c^4 - 60*a^5*b^6*c^3 - 340*a^6*b^2*c^6 + 364*a^6*b^3*c^5 + 136*a^6
*b^4*c^4 - 100*a^7*b^2*c^5))/c^8 - (((2048*tan(x/2)*(8*a*c^14 - 64*a^2*c^13 + 80*a^3*c^12 + 168*a^4*c^11 - 192
*a^5*c^10 - 136*a^6*c^9 + 72*a^7*c^8 - 2*b^2*c^13 + 6*b^3*c^12 - 17*b^4*c^11 + 33*b^5*c^10 - 49*b^6*c^9 + 61*b
^7*c^8 - 52*b^8*c^7 + 36*b^9*c^6 - 24*b^10*c^5 + 8*b^11*c^4 + 84*a*b^2*c^12 - 178*a*b^3*c^11 + 295*a*b^4*c^10
- 416*a*b^5*c^9 + 375*a*b^6*c^8 - 308*a*b^7*c^7 + 244*a*b^8*c^6 - 72*a*b^9*c^5 - 8*a*b^10*c^4 + 184*a^2*b*c^12
- 328*a^3*b*c^11 - 16*a^4*b*c^10 + 496*a^5*b*c^9 - 88*a^6*b*c^8 - 416*a^2*b^2*c^11 + 770*a^2*b^3*c^10 - 723*a
^2*b^4*c^9 + 779*a^2*b^5*c^8 - 732*a^2*b^6*c^7 + 80*a^2*b^7*c^6 + 112*a^2*b^8*c^5 - 8*a^2*b^9*c^4 + 180*a^3*b^
2*c^10 - 494*a^3*b^3*c^9 + 521*a^3*b^4*c^8 + 572*a^3*b^5*c^7 - 424*a^3*b^6*c^6 + 56*a^3*b^7*c^5 + 8*a^3*b^8*c^
4 + 234*a^4*b^2*c^9 - 1152*a^4*b^3*c^8 + 416*a^4*b^4*c^7 - 140*a^4*b^5*c^6 - 72*a^4*b^6*c^5 + 64*a^5*b^2*c^8 +
192*a^5*b^3*c^7 + 220*a^5*b^4*c^6 - 256*a^6*b^2*c^7 - 24*a*b*c^13))/c^8 - (((2048*(16*a^3*c^13 - 32*a^2*c^14
+ 176*a^4*c^12 + 176*a^5*c^11 + 48*a^6*c^10 - 2*b^4*c^12 + 6*b^5*c^11 - 18*b^6*c^10 + 26*b^7*c^9 - 12*b^8*c^8
+ 16*a*b^2*c^13 - 40*a*b^3*c^12 + 122*a*b^4*c^11 - 192*a*b^5*c^10 + 74*a*b^6*c^9 + 20*a*b^7*c^8 + 64*a^2*b*c^1
3 - 144*a^3*b*c^12 - 352*a^4*b*c^11 - 144*a^5*b*c^10 - 204*a^2*b^2*c^12 + 388*a^2*b^3*c^11 - 50*a^2*b^4*c^10 -
182*a^2*b^5*c^9 + 4*a^2*b^6*c^8 - 260*a^3*b^2*c^11 + 496*a^3*b^3*c^10 + 10*a^3*b^4*c^9 - 20*a^3*b^5*c^8 - 148
*a^4*b^2*c^10 + 116*a^4*b^3*c^9 + 8*a^4*b^4*c^8 - 44*a^5*b^2*c^9))/c^8 - (1024*tan(x/2)*(b^2*2i - a*c*2i + c^2
*1i)*(32*a*c^16 - 64*a^2*c^15 - 128*a^3*c^14 + 64*a^4*c^13 + 96*a^5*c^12 - 8*b^2*c^15 + 24*b^3*c^14 - 32*b^4*c
^13 + 32*b^5*c^12 - 24*b^6*c^11 + 8*b^7*c^10 + 144*a*b^2*c^14 - 200*a*b^3*c^13 + 184*a*b^4*c^12 - 56*a*b^5*c^1
1 - 8*a*b^6*c^10 + 288*a^2*b*c^14 + 352*a^3*b*c^13 - 32*a^4*b*c^12 - 320*a^2*b^2*c^13 + 8*a^2*b^3*c^12 + 96*a^
2*b^4*c^11 - 8*a^2*b^5*c^10 - 272*a^3*b^2*c^12 + 40*a^3*b^3*c^11 + 8*a^3*b^4*c^10 - 56*a^4*b^2*c^11 - 96*a*b*c
^15))/c^11)*(b^2*2i - a*c*2i + c^2*1i))/(2*c^3))*(b^2*2i - a*c*2i + c^2*1i))/(2*c^3))*(b^2*2i - a*c*2i + c^2*1
i))/(2*c^3))*(b^2*2i - a*c*2i + c^2*1i))/(2*c^3) + (((2048*tan(x/2)*(20*a*b^12 + 4*b^12*c - 4*b^13 - 40*a^2*b^
11 + 40*a^3*b^10 - 20*a^4*b^9 + 4*a^5*b^8 + 2*a^4*c^9 - 18*a^5*c^8 + 38*a^6*c^7 + 2*a^7*c^6 - 44*a^8*c^5 + 12*
a^9*c^4 + b^8*c^5 - b^9*c^4 + 4*b^10*c^3 - 4*b^11*c^2 - 8*a*b^6*c^6 + 4*a*b^7*c^5 - 31*a*b^8*c^4 + 20*a*b^9*c^
3 - 20*a*b^10*c^2 - 160*a^2*b^10*c + 320*a^3*b^9*c + 26*a^4*b*c^8 - 300*a^4*b^8*c - 84*a^5*b*c^7 + 136*a^5*b^7
*c + 2*a^6*b*c^6 - 24*a^6*b^6*c + 168*a^7*b*c^5 - 92*a^8*b*c^4 + 20*a^2*b^4*c^7 + 8*a^2*b^5*c^6 + 82*a^2*b^6*c
^5 + 6*a^2*b^7*c^4 + 8*a^2*b^8*c^3 - 44*a^2*b^9*c^2 - 16*a^3*b^2*c^8 - 40*a^3*b^3*c^7 - 104*a^3*b^4*c^6 - 132*
a^3*b^5*c^5 + 34*a^3*b^6*c^4 + 72*a^3*b^7*c^3 + 460*a^3*b^8*c^2 + 82*a^4*b^2*c^7 + 174*a^4*b^3*c^6 + 41*a^4*b^
4*c^5 - 149*a^4*b^5*c^4 - 660*a^4*b^6*c^3 - 900*a^4*b^7*c^2 - 90*a^5*b^2*c^6 + 96*a^5*b^3*c^5 + 541*a^5*b^4*c^
4 + 1156*a^5*b^5*c^3 + 764*a^5*b^6*c^2 - 204*a^6*b^2*c^5 - 704*a^6*b^3*c^4 - 840*a^6*b^4*c^3 - 300*a^6*b^5*c^2
+ 384*a^7*b^2*c^4 + 272*a^7*b^3*c^3 + 44*a^7*b^4*c^2 - 32*a^8*b^2*c^3 + 24*a*b^11*c))/c^8 + (((2048*(12*a^3*c
^11 - 28*a^4*c^10 - 44*a^5*c^9 + 72*a^6*c^8 + 88*a^7*c^7 + 12*a^8*c^6 + b^5*c^9 - 4*b^6*c^8 + 10*b^7*c^7 - 20*
b^8*c^6 + 29*b^9*c^5 - 30*b^10*c^4 + 26*b^11*c^3 - 12*b^12*c^2 - 6*a*b^3*c^10 + 27*a*b^4*c^9 - 72*a*b^5*c^8 +
154*a*b^6*c^7 - 238*a*b^7*c^6 + 251*a*b^8*c^5 - 228*a*b^9*c^4 + 98*a*b^10*c^3 + 20*a*b^11*c^2 + 8*a^2*b*c^11 -
68*a^3*b*c^10 + 112*a^4*b*c^9 + 100*a^5*b*c^8 - 200*a^6*b*c^7 - 96*a^7*b*c^6 - 47*a^2*b^2*c^10 + 145*a^2*b^3*
c^9 - 354*a^2*b^4*c^8 + 612*a^2*b^5*c^7 - 655*a^2*b^6*c^6 + 635*a^2*b^7*c^5 - 202*a^2*b^8*c^4 - 222*a^2*b^9*c^
3 + 4*a^2*b^10*c^2 + 239*a^3*b^2*c^9 - 524*a^3*b^3*c^8 + 536*a^3*b^4*c^7 - 564*a^3*b^5*c^6 - 115*a^3*b^6*c^5 +
856*a^3*b^7*c^4 + 2*a^3*b^8*c^3 - 20*a^3*b^9*c^2 - 37*a^4*b^2*c^8 + 9*a^4*b^3*c^7 + 583*a^4*b^4*c^6 - 1362*a^
4*b^5*c^5 - 152*a^4*b^6*c^4 + 156*a^4*b^7*c^3 + 8*a^4*b^8*c^2 - 399*a^5*b^2*c^7 + 904*a^5*b^3*c^6 + 394*a^5*b^
4*c^5 - 388*a^5*b^5*c^4 - 60*a^5*b^6*c^3 - 340*a^6*b^2*c^6 + 364*a^6*b^3*c^5 + 136*a^6*b^4*c^4 - 100*a^7*b^2*c
^5))/c^8 + (((2048*tan(x/2)*(8*a*c^14 - 64*a^2*c^13 + 80*a^3*c^12 + 168*a^4*c^11 - 192*a^5*c^10 - 136*a^6*c^9
+ 72*a^7*c^8 - 2*b^2*c^13 + 6*b^3*c^12 - 17*b^4*c^11 + 33*b^5*c^10 - 49*b^6*c^9 + 61*b^7*c^8 - 52*b^8*c^7 + 36
*b^9*c^6 - 24*b^10*c^5 + 8*b^11*c^4 + 84*a*b^2*c^12 - 178*a*b^3*c^11 + 295*a*b^4*c^10 - 416*a*b^5*c^9 + 375*a*
b^6*c^8 - 308*a*b^7*c^7 + 244*a*b^8*c^6 - 72*a*b^9*c^5 - 8*a*b^10*c^4 + 184*a^2*b*c^12 - 328*a^3*b*c^11 - 16*a
^4*b*c^10 + 496*a^5*b*c^9 - 88*a^6*b*c^8 - 416*a^2*b^2*c^11 + 770*a^2*b^3*c^10 - 723*a^2*b^4*c^9 + 779*a^2*b^5
*c^8 - 732*a^2*b^6*c^7 + 80*a^2*b^7*c^6 + 112*a^2*b^8*c^5 - 8*a^2*b^9*c^4 + 180*a^3*b^2*c^10 - 494*a^3*b^3*c^9
+ 521*a^3*b^4*c^8 + 572*a^3*b^5*c^7 - 424*a^3*b^6*c^6 + 56*a^3*b^7*c^5 + 8*a^3*b^8*c^4 + 234*a^4*b^2*c^9 - 11
52*a^4*b^3*c^8 + 416*a^4*b^4*c^7 - 140*a^4*b^5*c^6 - 72*a^4*b^6*c^5 + 64*a^5*b^2*c^8 + 192*a^5*b^3*c^7 + 220*a
^5*b^4*c^6 - 256*a^6*b^2*c^7 - 24*a*b*c^13))/c^8 + (((2048*(16*a^3*c^13 - 32*a^2*c^14 + 176*a^4*c^12 + 176*a^5
*c^11 + 48*a^6*c^10 - 2*b^4*c^12 + 6*b^5*c^11 - 18*b^6*c^10 + 26*b^7*c^9 - 12*b^8*c^8 + 16*a*b^2*c^13 - 40*a*b
^3*c^12 + 122*a*b^4*c^11 - 192*a*b^5*c^10 + 74*a*b^6*c^9 + 20*a*b^7*c^8 + 64*a^2*b*c^13 - 144*a^3*b*c^12 - 352
*a^4*b*c^11 - 144*a^5*b*c^10 - 204*a^2*b^2*c^12 + 388*a^2*b^3*c^11 - 50*a^2*b^4*c^10 - 182*a^2*b^5*c^9 + 4*a^2
*b^6*c^8 - 260*a^3*b^2*c^11 + 496*a^3*b^3*c^10 + 10*a^3*b^4*c^9 - 20*a^3*b^5*c^8 - 148*a^4*b^2*c^10 + 116*a^4*
b^3*c^9 + 8*a^4*b^4*c^8 - 44*a^5*b^2*c^9))/c^8 + (1024*tan(x/2)*(b^2*2i - a*c*2i + c^2*1i)*(32*a*c^16 - 64*a^2
*c^15 - 128*a^3*c^14 + 64*a^4*c^13 + 96*a^5*c^12 - 8*b^2*c^15 + 24*b^3*c^14 - 32*b^4*c^13 + 32*b^5*c^12 - 24*b
^6*c^11 + 8*b^7*c^10 + 144*a*b^2*c^14 - 200*a*b^3*c^13 + 184*a*b^4*c^12 - 56*a*b^5*c^11 - 8*a*b^6*c^10 + 288*a
^2*b*c^14 + 352*a^3*b*c^13 - 32*a^4*b*c^12 - 320*a^2*b^2*c^13 + 8*a^2*b^3*c^12 + 96*a^2*b^4*c^11 - 8*a^2*b^5*c
^10 - 272*a^3*b^2*c^12 + 40*a^3*b^3*c^11 + 8*a^3*b^4*c^10 - 56*a^4*b^2*c^11 - 96*a*b*c^15))/c^11)*(b^2*2i - a*
c*2i + c^2*1i))/(2*c^3))*(b^2*2i - a*c*2i + c^2*1i))/(2*c^3))*(b^2*2i - a*c*2i + c^2*1i))/(2*c^3))*(b^2*2i - a
*c*2i + c^2*1i))/(2*c^3)))*(b^2*2i - a*c*2i + c^2*1i)*1i)/c^3